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on  reverse  >t  necessary  and  identify  by  block  number) 

h  This  paper  concerns  techniques  for  processing  sonar  array  data  ij^nTer  to 
information  about  the  spatial  distribution  of  acoustic  sources.  The  goal1-'  .  f 
uresent.  a  method  for  determining  the  spectral  density  and  the  signal  wavefronts 
if  i.  sources  using  a  minimum  of  prior  assumptions.  Consideration  of  this 
t  is  the  next  logical  step  along  the  evolutionary  path  that  spatial  process¬ 
's  taken  over  the  last  few  decades.  This  path  is  characterized  by  a  decreas¬ 
'd.  of  confidence  in  any  kind  of  prior  knowledge  about  the  effects  of  the 
tiler.  medium,  which  has  proven  to  bo  more  and  more  /elusive  and  unpredictable 
set.  ion  ranges  and  array  size  have  increased.  •  . 
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(’J)  The  passive  sonar  problem  is  defined  and  the  traditional  approach  to  a 
solution  using  matched  filters  is  presented.  Maximization  of  the  directivity 
factor  for  known  signals  with  static,  nonadaptive  filters  emerges  as  the  prime 
goal  of  traditional  processing.  Care  is  then  taken  to  properly  introduce  the 
cross-spectral  density  matrix  which  plays  a  crucial  role  in  adaptive  processing. 

(y)  A  realization  of  adaptive  noise  cancellation  through  the  use  of  signal- 
free  references  is  achieved  with  the  development  of  a '^Correlo-Filtre" . 

Unfortunately  adaptive  noise  cancellation  leads  to  degradation  when  a  supposedly 
known  wavefront  deviates  from  an  expected  shape.  Robustness  and  the  need  for 
a  processor  which  does  not  require  prior  knowledge  of  the  wavefront  is  then 
discussed. 

fU)  A  novel  method  of  identifying  completely  'unknown  wavefronts  is  proposed. 
Using  the  eigenvalues  of  the  cross-spectral  density  matrix  allows  the  number  of 
sources  to  be  determined.  The  wavefront  shapes  are  then  shown  be  scrambled  within 
the  eigenvectors  by  an  unknown  unitary  matrix.  Additional  information  needed  to 
select,  a  specific  unitary  matrix  is  introduced  by  modeling  the  sources  as  a 
coherent  sum  oif  “conditional  plane  waves-.  While  the  existence  of  a  solution 
is  posel,  no  specific  solution  is  obtaineddbecause  of  the  large  computational  effort 
required.  \ 
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|)r.  Mmuuz  has  worked  in  the  area  of  acoustic  signal  processing  since  the 
md  of  World  War  II.  As  f lie  first  director  of  the  French  naval  laboratory 
(ironpe  Etudes  Rccherches  Detection  Sousmarines,  he  guided  the  develop¬ 
ment  of  French  acoustic  signal  processing  from  the  ashes  of  war  to  one  of 
ill*'  premier  laboratories  in  existence  today.  Now  retired,  this  paper  repre¬ 
sents  a  desire  of  Dr.  Mermoz  to  share  with  us  his  knowledge  of  the  past, 
understanding  of  the  present,  and  vision  of  the  future  in  passive  sonar  sig¬ 
nal  processing.  Hence  this  report  is  divided  into  three  main  chapters  (2-4) 
<••>  denied  with  each  evolutionary  stage  in  our  understanding  of  acoustic 
signal  processing. 

In  chapter  2  he  defines  the  passive  sonar  problem  and  presents  the 
traditional  approach  to  a  solution  using  matched  filters.  Maximization 
of  the  directivity  factor  for  known  signals  with  static,  nonadaptivc  filters 
emerges  as  the  prime  goal  of  traditional  processing.  He  then  takes  some 
rare  to  properly  introduce  the  cross-spectral  density  matrix  which  plays  a 
crucial  role  in  adaptive  processing  discussed  in  the  remainder  of  this  report  . 

In  chapter  3  he  turns  his  attention  to  contemporary  methods  for  the 
adaptive  cancelation  of  noise  through  the  use  of  signal-free  references.  He 
achieves  a  realization  of  noise  cancelation  with  the  development,  of  a  “Corre- 
|i>  Kiltie".  Unfortunately  adaptive  noise  cancelation  works  too  well,  leading 
t<>  degradation  when  a  supposedly  known  wavefront  deviates  from  what 
U  expected.  He  then  proceeds  to  discuss  robustness  and  the  need  for  a 
processor  which  does  not  require  prior  knowledge  of  the  wavefront. 

I  Inis  we  are  lead  into  chapter  4  in  which  Dr.  Mermoz  puts  forth  a  novel 
method  of  identifying  completely  unknown  wavefronts.  Using  the  eigenval¬ 
ues  of  the  cross-spect ral  density  matrix  allows  the  number  of  sources  to 
be  determined.  The  wavefront  shapes  are  then  shown  be  scrambled  within 
the  eigen vt<  l oi bv  an  unknown  unitarv  matrix.  Additional  information 


i 


needed  to  select  a  specific  unitary  matrix  is  introduced  by  modeling  the 
sources  as  a  coherent  sum  of  “conditional  plane  waves”.  While  the  exis- 
t<*irce‘-of  a  solution  is  posed,  no  specific  solution  is  obtained  because  of  the 


Targe1  Computational  effort  required.  This  then  is  his  vision  of  one  future 
area  of  profitable  research  in  sonar  array  processing. 

Sohie  mention  must  be  made  of  the  nonstandard  notation  used  to  dif¬ 
ferentiate  between  various  types  of  vectors  and  matrices.  Row  vectors  ac¬ 
cented  with  a  left  arrow,  e.g.  R,  are  used  throughout  this  paper  to  denote 
complex  filter  coefficients.  Column  vectors  accented  with  a  right  arrow, 
e.g.  S,  are  used  to  denote  a  complex  wavefront.  Matched  filters  tuned  to 
a  particular  wavefront  are  represented  by  the  complex  conjugate  transpose 

of  the  wavefront.  S ^ .  Matrices  are  accented  with  a  double  arrow,  C. 

Cross-spectral  density  matrices  (CSDM)  of  different  types  are  differ¬ 
entiated  using  subscripts.  Throughout  sections  2. 1-3.1  C  represents  noise 
alone,  while  in  the  remainder  of  the  paper  C  represents  a  transformed  noise 
alone  CSDM.  The  CSDM  ol  perfect  sources  is  represented  by  C,  while  the 
raw  data  matrix  is  denoted  by  Cr. 

The  Naval  Research  Laboratory  contracted  Dr.  Mermoz  to  write  this 
report  in  August  1983.  The  finished  report  was  delivered  to  NRL  on 
20  June  198(5.  Since  that  time  the  paper  has  been  revised  on  a  part-time 
basis  by  the  editor. 

Any  technical  correspondence  concerning  this  paper  should  be  addressed 
to  the  author  at:  982  Chemin  de  la  Calade,  83140  SIX-FOURS,  France. 
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Editor 


Acknowledgement 


M 


9 

I'll.-  author  expresses  his  thanks:  to  Dr.  Raymond  Fitzgerald  of  ONR 
who.  while  still  working  for  NRL.  suggested  this  eontract.  did  (In'  work 
9?  necessary  to  get  it  granted,  and  edited  chapter  2;  to  the  Office  of  Naval 

*.’•  Research  for  funding  this  work;  to  the  NRL  authoritit.es  who  agreed  to 

monitor  it,  particularly  Mr.  Burton  Hurdle;  most  of  all,  to  Dr.  Peter 
»■*  Mignerey  of  NRL  who  toiled  over  a  manuscript  in  approximate  English 

*“  and  whose  fruitful  criticisms  have  greatly  improved  the  presentation. 


Contents 


1  INTRODUCTION  1 

2  DIRECTIVITY  AND  SOURCE  IDENTIFICATION  5 

3  ADAPTIVE  SIGNAL  PROCESSING  WITH  A  KNOWN 

WAVEFRONT  19 

4  IDENTIFICATION  OF  SIGNALS  WITH  UNKNOWN 

WAVEFRONTS  41 

5  GENERAL  CONCLUSION  71 

A  Classical  Relation  between  Cross-Spectral  Densities  Before 
and  After  Filtering  77 

B  Matrix  of  Filters  and  Related  Equations  79 

C  Classical  Factorization  of  a  CSDM  -  Normalized  Eigen¬ 
vectors  83 

D  Orthogonal  Images  of  a  Set  of  Correlated  Noises  85 

E  Canonical  Form  of  a  CSDM  87 

F  Preliminary  "Theorem"  for  Section  4.0  89 

Bibliography  91 


List  of  Figures 


2.1  Simplified  Processor  with  (/V  -  1)  Decotrelafed  Signal  Free 

References .  27 

2.2  Processor  with  a  Single  SFR  Noise  Subtracter .  29 

2.2  Adaptive  Processor  with  a  Single  SFR  C'orrelo-Filter .  34 

3.4  C'orrelo-Filter  in  the  Time  Domain  (with  Sell  Whitening  of 

the  SFR  1 .  35 

3.5  General  Adaptive  Processor  with  (A’  -  1)  Signal-Free  Refer¬ 
ences .  36 

R.l  Matrix  of  Linear  Filters .  SO 


( 


Chapter  1 

INTRODUCTION 


This  paper  deals  with  teclnii<|ues  tor  processing  sonar  array  data  in  order  to 
obtain  information  about  the  spatial  distribution  of  the  acoustic  sources. 
I  he  goal  is  to  present  a  method  for  identifying  sources  whose  frequency 
spectrum  and  wavefront  at  the  array  are  unknown.  Consideration  of  this 
problem  is  the  next  logical  step  along  the  evolutionary  path  that  spatial 
processing  has  taken  over  the  last  few  decades.  Phis  path  is  characterized 
by  a  decreasing  level  of  confidence  in  any  kind  of  prior  knowledge  about  the 
effects  of  the  propagation  medium,  which  has  proven  to  be  more  and  more 
elusive  and  unpredictable  as  detection  ranges  and  array  size  have  increased. 

In  order  to  provide  a  proper  introduction  for  this  problem,  it  is  useful  to 
retrace  the  evolutionary  path  with  a  modern  point  of  view  so  that  the  new 
concepts  can  be  progressively  introduced.  As  a  result,  this  introduction  is 
divided  into  three  sections,  which  correspond  roughly  to  the  three  major 
steps  in  the  evolution  of  spatial  signal  processing. 

The  concept  of  directivity  index  has  dominated  (if  not  choked)  spatial 
processing  lor  decades.  In  fact,  the  ut  ility  of  direct  ivity  index  as  a  descriptor 
requires  very  stringent  conditions  to  be  satisfied:  both  the  signal  wavefront 
and  I  lie  spatial  structure  of  the  noise  field  must  be  known.  The  signal 
wavefront  is  usually  assumed  to  be  a  plane  wave,  which  implies  that  (he 
propagation  medium  is  homogeneous.  I  lie  noise  field  is  usually  assumed 
to  be  "oiunidirec  I  ional.  which  again  implies  a  homogeneous  medium  and. 
in  addition,  completely  determines  the  noise  correlation  as  a  function  of 
the  spatial  separation  between  array  sensors.  Medium  behavior  is  therefore 
frozen  in  the  notion  of  directivity  index:  the  only  option  loll  for  spatial 
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processing  is  to  try  to  improve  the  signal  to  noise  ratio,  which  depends  on 
the  number  ol  sensors  and  on  their  geometry. 

Directivity  index  was  a  satisfactory  descriptor  for  the  high-frequency 
short-range  sonars  of  an  earlier  generation.  Many  senior  scientists,  however, 
can  remember  a  time  when  the  sonar  equation — in  decibels  was  allowed 
an  additional  term,  A.  for  "anomaly."  Because  the  medium  was  frozen  and 
allowed  tor  no  abnormalities,  the  concept  had  to  evolve. 

Circa  !!><>()  it  became  obvious  I  hat  a  description  ol  t  he  noise  field  as  om¬ 
nidirectional  was  no  longer  satisfactory.  Moreover.it  became  apparent  that 
I  he  description  of  t  he  spatial  struct  ure  as  frozen  was  to  be  questioned.  The 
new  picture  ol  a  time- varying  spatial  structure  called  for  the  development 
o!  t ime- varying  processors.  This  led  to  the  introduction  of  adaptive  spatial 
processing  and  adaptive  beam  forming,  in  which  a  n //  thought  of  usimj  prior 
kimtrli  iltjt  uhoiit  thi  noist  Ji<  Id  if  dismissed.  Instead  the  description  of  the 
noise  field  reflected  in  the  time-dependent-  correlation  matrix  changes  with 
time.  Here  the  notion  of  correlation  matrix  appears  explicitly,  whereas  it 
was  only  implicit  m  the  concept  of  directivity  index.  Although  the  corre¬ 
lation  matrix  plays  an  important  role  in  adaptive  beamforming,  that,  role 
remains  limited  compared  to  what  it  will  be  in  the  approach  presented  in 
chapter  i. 

Although  the  notion  ol  a  predetermined,  staf ic  noise  field  had  been  re¬ 
laxed,  the  signal  wavefront  was  presumed  still  to  be  a  plane  wave  whose 
direction  ol  propagation  was  known.  The  resulting  beam  former  operated 
to  look  simultaneously  in  all  possible  directions  with  preformed  beams.  In 
chapter  •>  we  will  generalize  (he  problem  slightly  by  assuming  that,  al¬ 
though  the  signal  wavefront  is  known,  it  is  not  necessarily  a  plane  wave. 
"Beamforming"  then  becomes  adaptive  processing  matched  to  a  given  sig¬ 
nal  wave-front. 

Adaptive  beamforming  was  doubtless  a  significant,  step  forward  in  that 
it  eliminated  all  prior  assumptions  about  the  structure  of  the  noise  held. 
Its  cost  lies  in  the  need  for  a  significant  amount  of  real  time  computing 
power  and  the  techniques  for  its  implementation  are  still  developing. 

It  was  predictable,  and  it  has  already  been  shown,  that  adaptive  beam- 
lorming  even  in  the  generalized  form  described  in  chapter  .!  would  be 
limited  by  the  only  assumption  it  retained:  namely,  that  the  signal  wave- 
front  is  known.  At  first  glance,  any  attempt  to  eliminate  this  last  assump¬ 
tion  would  seem  hopeless.  Nevertheless,  a  deeper  analysis  ol  the  structure 
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of  i  i  t  <  -  i.  .it  flat  i<  >ii  matrix.  shows  that,  iillhotn/h  .••nun  prior  mforniu  I  ion  is 
nijinrtti.  tli('  amount  is  less  than  that  contained  in  a  complete  description  of 
tin*  signal  wavefront.  Moreover,  assumption  ol  ;i  too  detailed  prior  descrip- 
i  i  i . :  -  oi  the  signal  wavetronl  can  easiK  prove  t  <  >  lie  incompatible  with  the 
measured  corn-  .it ion  matrix.  Ihe  temptation  alwaxs  exists,  in  that  case, 
to  hliime  tin  i  .  i  an  pat  il  nht  \  on  inaccnrac  ics  in  I  lie  <>l  iniiil  n  I  correlation 
main.'..  I  ha’  tins  is  not  the  source  o!  the  problem  is  suggested  both  hy 
("ir  inci'easiin-'  abilitv  to  ar<|'iire  accurate  data  and  b\  i  he  do  elopuient  *  >  { 
i  (‘lined  est  i:  •  i  a  i  .i .  i .  methods.  We  shall  comment  onl\  t  •  r  i « ■  i  I  \  on  the  eopi 
oils  lit  era  i  n  re  •!(,'oted  to  estimation  methods:  essentialh  we  slndl  accept 
mid  1 1  ust  the  methods  its  they  exist.  Having  done  so  we  can  consider 
at  temp!  inu  to  i  opr  with  the  complexity  of  the  medium  and  its  effects  in 
older  to  determine  hot  i  the  signal  wavefront  and  its  spectral  density  with 
I  a  rue  iirne.s.  |  hi-  is  still  adaptive  processing,  but  not  simply  benmlorming. 
and  it  require-.  a  huge  amount  of  real  time  computing  power. 

What  we  are  attempting,  to  do.  then,  is  to  present  a  method  lor  the 
"nit  ntilirnhoo"  li.e.  determination  ol  the  spectral  density  and  the  signal 
■a  a  ve  Iron  I  si  ol  one  or  more  sources  using  a  minimum  of  prior  assn  nipt  ions. 
I  In  \ ( ■  r \  notions  of  siu.ua!  and  noise  will  be  slightly  dillevent  from  what 
’  lie  \  wel  l  ’  I  :  1  >  i:.l  pi  e|  s  2  a  ltd  -i. 

We  shad  not  -leal  with  the  problem  ol  source  localization  nor  shall  we 
,|.;d  with  am,  -p- i  ihe  model  of  I  he  medium.  Instead  we  shall  indicate  in 
chapter  I  v.li.il  can  be  done  in  approaching  the  problem  of  source  loealiza 
tii, n  (i.e.  del.  i  mmim:  the  coordinates  of  the  source).  In  previous  papers, 
most  of  which  tii"  in  l-'remh.  we  mixed  the  problems  of  identification  and 
|o,  ali/alion  i  -  i.  u-millx  done  in  more  supeifieial  treatments.  What  we 
dndl  siio i'  here  i-  tiiat  localization,  and  only  localization,  is  inseparable 
limn  soiiu  sort  of  description  or  model  of  the  medium.  Separating  the 
i  wo  probh  no  .'.ill  slu'd  clarity  on  both.  W  hile  the  problem  of  wavefront 
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computing  power  before  any  real-time  implementation  becomes  possible. 
Nevertheless,  some  experimental  investigation  is  possible  today  and  those 
particular  aspects  ol  the  theory  which  require  early  verification  will  be 
indicated.  Indeed  the  kind  of  "weak"  assumptions  which  remain  in  chapter 
1  are.  in  fact,  the  last  ramparts  in  spatial  processing.  II  these  01  some 
similar  assumptions  were  totally  unacceptable  (or  a  given  problem,  the 
implication  would  be  that  no  sensor  combination  would  be  more  ellective 
t  han  a  single  sensor.  Such  a  eatast rophic  case  seems  quite  unlikely,  however, 
because  arrays  have  t  hits  far  proven  to  be  useful  tools  which  do  provide  some 
gain,  although  generally  not  the  gain  predicted  by  current  theories. 
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DIRECTIVITY  AND 
SOURCE  IDENTIFICATION 


2.1 


Much  lias  been  written  about  directivity,  yet  the  concept  actually  contains 
t  wo  different  not  ions.  The  more  familiar  not  ion  of  angular  resolving  power, 
which  is  related  to  beamwidth,  has  significance  for  both  transmitting  and 
receiving  arrays.  The  scientific  literature  dealing  with  the  minimization  of 
beamwidth  and  the  reduction  of  sidclobe  levels  is  extensive.  Imr  receiving 
arrays  it  is  well  known  that,  by  using  nonlinear  combinations  of  sensor  out¬ 
puts,  one  can  obtain  narrower  main  beams  and  hence  higher  resolution  t ban 
with  linear  techniques,  at  least  for  strong  signals.  Such  nonlinear  processing 
techniques  are  inetlicient  and  unstable  and  can  lead  to  erroneous  results  for 
weak  signals  in  a  noisy  environment.  The  second  notion  contained  in  the 
concept  of  directivity  is  that  of  directivity  factor  or,  in  decibels,  directivity 
index,  which  is  a  measure  of  the  spatial  array  gain  (i.e.  the  increase  in  the 
signal -to-noise  ratio)  for  signals,  weak  or  strong,  in  the  presence  of  noise. 

Because  our  intention  is  to  treat  the  detection  and  identification  of  weak 
signals  in  background  noise,  we  are  concerned  here  only  with  the  second 
notion  o!  directivity  and  with  linear  combinations  <>l  sensors.  Tims,  we 
deal  with  processors  in  which  each  sensor  output  is  passed  through  a  linear 
filler.  I'.acli  filter  can  be  described  in  terms  either  <>l  its  impulse  response 
or  ol  its  complex  frequency  response,  which  form  a  homier  pair.  The  term 
"combination'  will  be  used  in  this  paper  to  refer  to  the  sum  of  the  outputs 
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of  these  linear  filters. 

In  a  nonadaptive  heamlornier.  such  as  one  having  classical  preformed 
beams,  the  filters  do  not  change  with  time;  in  adaptive  beamlormers  I  lu- 
filter  characteristics  are  continually  adjusted  in  accordance  with  the  in  lor 
illation  obtained  from  estimates  of  the  evolving  noise  correlation  matrix. 
The  permanence  of  the  filters  in  classical  beamforming  is  a  direct  rouse 
quence  of  tin1  very  rigid  assumptions  made  about  the  structure  of  (  he  signal 
and  noise.  These  assumptions  are: 

1.  The  signal  is  uncorrelated  with  the  noise  and  is  carried  by  a  single 
plane  wave  coming  from  a  direction  specified  by  ^  which  stands  lor 
the  two  angles,  bearing  and  tilt. 

2.  The  noise  is  omnidirectional  (i.e.  spherically  isotropic).  This  means 
that  for  each  frequency  the  noise  field  in  the  vicinity  of  the  array 
is  composed  of  a  superposition  of  uncorrelated  plane  waves  which 
can  be  described  as  carrying  equal  power  from  all  directions  or  as 
being  uniformly  distributed  over  all  solid  angles  du’.  Implicit  in  both 
assumptions  is  the  notion  that  the  medium  is  isotropic. 


For  the  signal  and  noise  fields  described  above,  the  behavior  of  the  array 
is  completely  described  by  its  plane  wave  response  function.  For  a  given 
sensor  combination  t  he  plane  wave  response  is  the  complex  function  /)( u.\  / ) 
which  gives  the  amplitude  and  phase  at  the  beamformer  out  put .  when  the 
input  is  a  unit  amplitude  plane  wave  having  frequency  /  and  coming  from 
the  direction  u>.  The  output  phase  is  measured  relative  to  one  of  the  sensors 
or  relative  to  some  reference  point  in  the  array. 

For  a  nonrandom  signal  coming  Irom  the  direction  u?o,  the  system  re 
spouse  at  frequency  /  is  ft{f)l)(usn, /).  where  <t(f)  is  the  complex  amplitude 
of  the  Fourier  component  at  frequency  /  measured  at  the  signal  origin.  For 
this  signal  the  power  spectral  density  at  the  output  is  s( /)!“  /)( u.-n .  / I " 
where  the  vertical  bars  denote  the  modulus  of  a  complex  quantity,  so  that 
l/J|2  DD'  with  the  star  denoting  complex  conjugation. 

II  the  signal  coming  Irom  the  u.-,,  direction  is  stochastic  with  powei  spec 
t  ral  density  »/(/).  then  the  power  sped  ml  density  at  the  beam  lormer  out  pul 
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where  l\  is  the  number  of  sensors.  This  set  of  numbers  is  the  transla 
lion  through  the  array  of  the  incoming  plane  wave  and  can  be  consid¬ 
ered  to  Ire  a  vector  in  an  A’-dimensional  complex  space.  For  any  given 
array  geometry  this  vector  is  easily  constructed.  Let  /»„(«.'./)  be  the 
Fourier-transformed  output  of  the  nth  sensor  evaluated  at  frequency 
/  and  let  //(u.\/)  be  the  vector  whose  A  components  are  the  h„. 

2.  The  second  stage  of  processing  is  performed  by  the  combination 
that  is,  by  summing  the  outputs  of  the  A  linear  filters.  These  filters 
are  described  by  their  frequency  responses.  )•„(/ ),  which  are  the  com¬ 
ponents  of  a  vector  R.(/). 

The  beamformer  out  put  resulting  from  these  two  stages  of  processing  is 

,Y 

DM/)  ---  V  f-.M 

I)  :  1 

Understanding  that  all  quantities  are  functions  of  frequency,  we  shall  hence 
forth  suppress  the  frequency  variable  in  order  to  simplify  the  notation. 

If  R  is  taken  to  be  a  row  vector  and  H(u’)  a  column  vector,  then  (2.1  ) 
becomes  the  scalar  product 

DM  -  RHM  (2.2) 

The  column  vector  HM  is  a  j)articular  case  of  what  we  shall  call  a  source 
vector"  corresponding  to  a  remote  source  transmitting  white  noise  which 
arrives  at  the  array  from  the  direction  u?.  The  corresponding  normalized 

vector  is  F(iv)  HM/|HM|  where  (HMl  -  |H*(sv  )H(u>  )|‘ /2  and  the 
superscript  ^  indicates  the  conjugate  transpose.  We  shall  refer  to  the  vector 
F  as  the  “wavefront  vector,”  although  if  is  not  the  physical  wavefront  (here 
a  plane  wave)  but  its  translation  through  the  array  of  sensors.  In  the 
concept  of  directivity,  it  is  implicitly  assumed  that  this  wavefront  vector  is 
perfectly  known  for  plane  waves  arriving  at  the  given  array  geometry. 

The  filter  response  r„  must  be  chosen  to  maximize  (he  ratio 
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Optimization  of  />( u.'o)  usually  leads  to  a  main  beam  pointed  in  the 
direction.  II.  however,  additional  constraints  are  placed  on  the  solution. 
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elicit  as  equality  ol  side-lobe  levels  or  a  predetermined  rale  of  decrease  in 
side  lobe  level,  then  the  maximum  response  of  the  array  may  not  be  in 
the  direction  ol  the  source.  In  other  words,  optimization  of  the  directivity 
factor  does  not  necessarily  meet  the  requirements  suggested  by  a  desire  to 
improve  angular  resolution.  Still,  optimization  of  the  directivity  factor  does 
provide  the  maximum  signal-to-noise  ratio  under  the  given  assumptions. 
This  is  one  illustration  ol  the  previously  mentioned  duality  in  the  concept 
o|  direct  ivil  v. 


W  ith  />(«.')  given  by  (2.1).  the  integral  in  the  denominator  of  p(u.’0)  can  be 
written  as 
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which  can  be  rewritten  as 


/  i  /l(a’  )|"  fru.' 

■hi 


,V  .V 


1  ’1 1 III  /  h  >1 

■hi 


(tt’  )/i„,(<^ )  <*’<- 


T  lie  t|ua  lit  it  y  )//  (u.1 ).  which  is  the  element  in  the  n"1  row  and  ni"' 

column  of  I  he  elementary  matrix  H(  u>  )Hf(  u.' ),  is  the  cross-spectral  density 
of  the  held  coming  from  the  directions  u.’.  This  can  be  seen  by  recalling  the 
general  rule: 

II  a  source  with  spectral  density  rr  passes  through  two  filters 
h„  and  h,„  in  parallel,  then  the  cross-spectral  density  of  the  two 
outputs  is  I],,,,,  n h „ li m  (see  appendix  A). 

In  the  particular  case  under  consideration  rr  1.  because  H  corresponds 
to  a  unit  amplitude  plane  wave.  Thus,  the  cross-sped ral  density  matrix 

for  a  point  like  source  in  the  far  field  is  a  dyad  of  the  form  HI  a.-  )H  1  ( u.’ ). 
In  older  to  understand  the  approach  in  tin'  later  sections  ol  this  paper,  it 
is  essential  to  observe  that  the  rank  of  this  matrix  is  1.  even  though  its 
order  is  A  .  Indeed,  every  vector  X  that  is  orthogonal  to  H  in  an  (A  I  )- 

dimensional  subspace  is  transformed  by  HH^  to  a  null  vector.  Obviously, 
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with  H  and  X  as  column  vectors,  we  have  H'X  -  0  which  implies  that 
(HH^)X  =  0.  The  only  nonzero  eigenvalue  of  HH  ^  is  the  scalar  H^H  and 
the  corresponding  eigenvector  is  H  itself. 


Looking  now  at  the  integral  on  the  right  hand  side  of  (2.0).  we  write 


Jn 


This  is  clearly  the  (mm)'1'  element  in  the  sum  of  an  infinite  number  of 
elementary  cross-spectral  density  matrices  linked  to  uncorrelated  point -like 
sources.  Therefore,  this  sum  is  the  (»,»?)"'  element  of  the  cross-spectral 
density  matrix  C  for  the  normalized  omnidirectional  noise.  Although  each 
elementary  contribution  to  matrix  C  is  a  rank-one  matrix,  the  rank  of 
C  is  not  unity.  Neither  is  it  less  than  the  order  N,  because  the  rank 
increases  every  time  an  elementary  contribution  is  added,  and  there  is  an 
infinity  of  such  contributions  (tins  will  be  better  illustrated  in  chapter  I ). 
As  a  result,  the  rank  and  order  of  C  are  A’  and  C  is  invertible  (i.e.  C  ’ 
exists).  Something  similar  would  happen  if  the  sensor  outputs  consisted 
of  independent  noise  with  the  same  spectral  density.  The  cross-spectral 
density  matrix  would  again  be  of  rank  A7,  but  it  would  have  a  simpler 
form,  being  proportional  to  a  unit  matrix  (which  is  not  line  for  matrix  C 


It  is  obvious  from  (2.0)  that  c„m  cm„.  The  matrix  C  is  llermitian, 

C  =  C^,  and  as  a  result  has  real,  positive  eigenvalues,  as  does  every  cross- 
spectral  density  matrix.  The  same  is  true  for  the  inverse  matrix  C  The 
right  hand  side  of  (2.-ri)  can  now  be  written  as 

_ I  \  * 

RCR  1  \  \  e„„, 


which  is  the  quadratic  form  for  the  row  vector  R  through  the  matrix  O. 
This  real,  positive  scalar  is  t  lie  sped  ral  density  "I  I  he  omnidirect ional  noise 
at  the  output  ol  the  combination  ol  sensors. 
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Looking  now  to  the  numerator  of  />(u>0)  in  (2.3), 

|/J(u.’o)('  —  D(u.'a)L)  (u>o) 


we  ran  use 


to  write 


r>( 
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V  h„("-'o)r„, 


N  AT 


D{a.'o)D  (a’o )  —  }*  ^  r?ir,„/bi(u-’o)/fm(u'’o). 

=  1  m  —  1 

The  quantity  /?„(u-’o  )/)m(u,'0)  is  (lie  element  in  the  n,fl  row  and  nitU  column 
ol  the  matrix  HnH(|  related  to  the  incoming  plane  wave  signal.  The  rank  of 

this  matrix  is  unity.  The  only  nonzero  eigenvalue  is  the  scalar  H,;H0  and 
tin-  associated  eigenvector  is  H0-  As  we  shall  see  later.  H0  is  (lie  source 
vector  of  a  remote  source  transmitting  a  white  noise  from  the  direction  u.>o. 
We  have 


|/dU>)|2  -  RHohJr1  -  |RH„ 


(2.7) 


which  is  simultaneously  the  quadratic  form  of  the  vector  R.  through  the 

matrix  HnHn  anti  the  stpiare  of  the  modulus  of  the  scalar  product  RH0. 

Maximizing  p(u.’n)  in  (1.3)  then  consists  ol  choosing  the  filters  R  so  as 
to  maximize 


p(uJo)  — 


|RHn| 


(2.S) 


RCR^ 

This  is  exactly  the  same  problem  that  is  encountered  in  determining  the 
spatial  matched  filler  for  a  nonrandom  signal  in  the  presence  of  noise,  whose 
complex  spectral  components  at  the  sensor  outputs  are  H(l.  with  a  cross- 
spectral  density  matrix  given  by  C. 

The  solution  is  well  known  and  can  he  obtained  in  several  ways.  One  of 
them  is  briefly  outlined  here.  In  order  to  present  it  in  classical  form,  let  us 

change  notation  by  replacing  the  row  matrix  R  with  the  row  matrix  X ^ . 
Ihe  column  matrix  X  then  represents  the  filters  through  I  lie  convention 
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x„  =  r'n .  In  the  new  notation  (2.8)  becomes: 


p{  ) 


XfH0 


X'CX 


Matrix  C,  being  Hermitian  and  invertible,  can  be  written  in  factored  form 
as  C  =  M^M  with  M  an  invertible  mat  rix.  Hence: 

C-'  =  M”1  ( K/rt \  -  lvf1  t' lvr ~ 1 1  ^  o  mi 


"‘(m')  =  1VT1  (iVT1) 


Defining  the  column  vectors  Y  and  P0  by: 

Y  =  MX, 

Po  =  (mM 


equation  (2.9)  becomes: 

-t-  2 

Y'Po 

p{uj0)  -  zirz  • 

Y’Y 

Now  the  Schwartz  inequality  states  that  for  any  Y  and  P0 

lvtp0‘  <  (yty)  (pjp0)  . 


Equality  occurs  when  Y  and  P0  are  parallel: 

Y  =  oP0.  (2.11 ) 

with  o  an  arbitrary  frequency  dependent  complex  scalar.  From  (2.12)  and 
(2.13),  the  maximization  of  p(u>0)  is  reached  through  (2.14).  This  maximum 
in  conjunction  with  (2.10)  and  (2.11  )  allow  X  and  R  to  be  written  as: 


aC'  H0, 
o  "h|c  1 : 


(2.15) 


which  is  a  spatial  generalization  of  the  matclied  filter  theorem.  The  max- 
imimum  value  of  p(u;(1)  is  then. 


PoPn 

H,tc  'H„. 


(2.10) 
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l  ln'  lilte-rs  R  ami  X  defined  by  (2.15)  are  multiplied  by  an  arbitrary 
lilte-r  factor  which  docs  not  alter  (provided  the  filter  has  no  null  in 

ils  lrci|ueney  response  across  the  useful  bandwidth).  11  some  of  the  filters 
arc  noucansa!  it  is  possible  lo  make  them  so  by  finding;  a  time  delay  (or  a 
phase-  filter  in  series:  exp|2/r/r] ).  which  translates  all  the  impulse'  responses 
to  causality.  Morcovcr.it  is  also  possible  to  get  an  inuiliord  signal  speed ral 
density  on  the  output  of  the  filter  combination.  Since  the  ultimate'  goal  is 
to  analyze  the  spectral  density,  an  unalteree!  signal  is  gene-rally  reepiiivel. 
( )nc  has  to  choose'  iRuHn!  -  1.  w  Inch  is  ecpiivalent  to 

|  / 1  ( u-’o .  / )  |  I ,  I  2. 1  7  ) 


because  eil  (2.7).  ( 'onse-epieuit  !v. 

In!  -  (rijc  ’Ho) 

1  pmtl.r  • 


(2. IS) 


A  little-  more  restrictive  condition  D(u>o,f)  -  1,  lor  all  p,  generally  makes 
the-  metre  familiar  directivity  functions  all-pass  filters  for  the  signal.  Any¬ 
way  the-  final  expressiem  for  X  and  R  after  (2.15)  and  (2. IS)  is,  within  an 
m  bit  rary  phase-  factor: 


X 

R 


C  !H(, 

hJc  ’h 

hJc  1 

H  l  C  II 


(2.19) 


2.7 

In  this  se-ction  we  are  cenicerned  with  directivity  and  omnidire-c  t ional  noise-. 
Vector  H  is  t  lie  array  respenise  to  a  unit  amplit  ude-  plane-  wave-  ( t  lie-  ge-ome-t  ry 
ol  tin-  array),  while-  C  eh-notes  a  uormaliz<-d  cross-spe-ct  ral  ele-nsity  matrix 
linked  to  H  by  (2.5).  It  is  note-worthy  that  all  relations  from  (2.9)  to  (2.19) 
are-  v  a  lid  for  any  v«-et  e  .r  H  a  ami  any  cre>ss-spect  ral  elensity  matrix  C.  for 
both  a  nonrandom  impulse-  signal.  »*>(/)  “carrie'il  by  the  wave-front'  Ho.  and 
a  st oe  liast  ie-  while-  signal  with  a  unit  spec!  ral  eh-nsit y  alsei  e  a r ri«-*l  by  H„,  the- 
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optimal  sot  o!  fillers  is  given  l»y  ( 12 . 1  f ) )  when  the  noise  spatial  slimline  i- 
representeil  by  the  non  normalized  matrix  C.  It  happens  that  w<  have  also 
parllv  treated  the  processing  problem  for  any  spatial  signal  represented  In 
a  given  vector  in  a  noise  structure  represented  l>y  a  given  matrix.  I  his  will 
I >e  used  in  chapt er  >■ 

Coming  back  to  directivity,  we  now  have  an  output  with  the  unaltered 
signal  sped  ral  density  ,/  mixed  with  the  fraction  of  t  he  <>m  indirect  ion  a  I  noise 
which  is  let  through  by  the  combination  R.  or  X.  According  to  section  Lhd. 
this  noise  spectral  density  is: 


h  /  CU’bdC 

mi 


/»( R  CR 


MX^CX) 


/>(  H„C  H(l ) 


I  lie  factor  Ho  C  'H()  is  unimportant  because  it  is  considered  to  be  perfectly 


known  in  lh<'  rlirrrl  i\ily  concept.  Since  the  signal  and  noise  are  umone 
lated.  the  total  spectral  density  on  the  output  is  ,/  j  h/  I  he  error 

in  I  he  spectral  analysis  of  the  signal  is  h  and  may  be  negligible  lor 

a  large  value  of  .  i.e..  a  large  array  gain.  Doing  belter  requires  some 

knowledge  of  l>  in  order  to  ’‘identify"  <1.  the  only  parameter  which  remains 
unident ilied  in  t lie  direct i vity  concept. 


file  real  problem  is  to  get  an  estimate  for  />.  Herein  arises  the  notion 
of  "signal- free  reference."  The  possibility  ol  building  a  combination  <>l 
sensors  from  which  only  the  noise  would  appear  is.  at  first  sight,  relative!;, 
easv.  If  we  choose  a  set  of  filters  represented  by  a  row  vector  V  such  that 
VII, ,  0.  linn  iiie  signal  with  a  wavefront  II,,  will  not  appear  at  lie 
output.  I  he  oi  ni  i  idi  red  ioua  1  noise  will  appear  at  the  beamloine'r  output 

with  spectral  densitv  h  V(A'  1  )  which  can  be  measured  and  Item  v  Inch 
/,  can  be  determined  In  principle  the  choice  o|  V  is  arbitrarv  v.ilhin  tin 
l  A  It  dimensional  -ubspace  orthogonal  to  II,,. 
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In  tin-  current  classical  preformed-beam  techniques,  (lie  role  ol  the 
"signal -lire"  combination  is  approximately  played  by  the  mean  value  of  the 
nut  puts  on  beams  adjacent  to  or  close  by  the  one  being  tested.  The  lact 
that  only  neighboring  beams  are  considered  clearly  reveals  our  uneasiness 
about  the  omnidirectionality  assumption. 

The  drawback  with  the  previous  procedure  is  that  there  may  be  some¬ 
thing  other  than  the  signal  o.*0  and  the  omnidirectional  noise  present.  In 
particular,  there  may  be  other  sources  which  are  not  well  described  in  the 
directivity  concept  as  unwanted  omnidirectional  noise.  These  other  sources 
can  add  unknown  contributions  to  both  combinations  R0  and  V.  and  can 
be  labeled  as  "unwanted  signals. 

Concerning  the  presence  of  unwanted  signals  arriving  at  the  array  as 
plane  waves  which  are  uncorrelated  with  each  other,  with  the  u.\,  source, 
and  with  the  omnidirectional  noise,  two  situations  can  arise: 

I.  In  the  first  situation,  the  direction  of  each  unwanted  signal  is  known, 
so  that  we  have  knowledge  of 

u.* |  and  therefore  IT 

and  therefore  HU 

u.*/'  and  therefore  H/-. 

Oik*  must  then  seek  a  “source- free*'  reference  for  t  In*  omnidirectional 
noise  which  has  a  vector  V  orthogonal  to  all  the  vectors 

Ho.H, . H, . H, 

for  this  to  be  possible,  there  must  exist  a  subspace  “free  of  sources.*’ 
I  1 1 1 1 s .  .V.  the  number  of  sensors  in  the  array,  must  be  larger  than 
/'  •  1 .  1  lie  number  of  sources  present .  This  is  but  t  he  first  illust  ration 
of  how  the  arrav  has  to  be  large  enough  to  surmount  the  complexify 
of  the  source  system.  This  topic  will  be  discussed  further  in  later 
sect  ions. 

||  the  condition  A  /’  •  1  is  fulfilled,  it  is  possdih*  to  remove  the 
omnidirectional  noise  spectral  density  from  tin*  outputs  ol  the  com¬ 
binations  deri  veil  I  rom  tin*  1I(.  according  to  I  2. 1  It.  \  t  the  H,,  output 
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we  will  still  have  spectral  contributions  due  to  the  unwanted  signals. 
Spectral  analysis  will  yield 

,t0  —  d  f  (I i  j R-o H l ! 2  +  •••  f  dp  |RoHp|2  at  the  H0  output, 

.4  i  -  d  |RiH0i‘  +  d\  +  •  •  •  +  dp  |RiHp|2  at  the  Hi  output, 

,4p  -  d  |RpHo|“  +  •  ■  •  +  dp  at  the  Hp  output. 

From  t  his  system  of  (  P  +  l )  equal  ions  for  the  ( P+  1 )  unknown  spectral 
densities  one  can  identify  all  the  sources  in  terms  of  their  individual 
spectral  densities. 

2.  In  the  second  situation,  the  directions  of  the  unwanted  sources  are 
unknown,  but  the  number  of  sources  is  expected  to  be  much  less  than 
.V.  It  is  then  still  possible,  in  principle,  to  remove  the  omnidirectional 
noise  from  the  H0  output.  By  sweeping  with  a  variable  combination 
V  all  the  subspace  orthogonal  to  H0.  we  can  determine  a  cone  in 
which  the  noise  is  both  minimum  and  constant  at  the  output  of  V. 
'{'his  value  is  qualified  to  represent  the  omnidirectional  noise.  Such  a 
process  is  certainly  painful  and  will  be  used  reluctantly,  because  the 
strict  omnidirectional  noise  assumption  can  not  be  trusted  to  such  an 
extent.  Anyway  there  is  now  no  way  to  eliminate  the  contributions 
from  a  random  distribution  of  other  sources.  Identification  of  the 
spectral  density  <1  is  hazardous  if  not  impossible. 

In  the  practical  implementation  of  classical  beamforming,  where  every 
wavefront  is  taken  to  be  a  plane  wave,  one  seeks  to  build  the  array  large 
enough  so  that  with  its  many  narrow  beams  it  can  minimize  the  scalars 
(R;,H,  )  and  limit  the  interference  between  sources.  Within  the  limited 
framework  provided  by  the  directivity  concept,  there  is  little  more  that  can 
be  done. 


2.9 

The  discussion  of  the  directivity  concept  can  be  summarized  as  follows: 

1.  Onlv  its  expression  as  a  receiving  array  train  need  be  considered  here. 


nun<  -min 


. SA 


ft  * 


Tin'  di'sc ■  r i | >1  ion  it  employs  of  signal  and  noise  st  rnel  me  is  over-simpli¬ 
fied  compared  to  that  in  the  real  world.  Thus,  this  description  relies 
on  extremely  stringent  assumptions. 


|  lie  st  at  us  of  waves  differing  from  t  he  one  carrying  t  lie  signal  is  some¬ 
what  vague  and  the  possible  presence  ol  unwanted  signals  consider- 
abl\  complicates  the  identification  of  the  wanted  signal. 


I .  I  lie  role  pin  veil  bv  t  he  correla  t  ion  matrix  is  modest .  Actually,  because 
the  function  />(.<.•./ !  is  usually  directly  accessible,  the  appearance  ol 
t  lie  cross- sped  ral  density  mat  rix  in  the  present  at  ion  in  sect  ion  2.1  was 
uunecessarv  and  somewhat  forced.  Die  only  reason  lor  introducing 
it  was  to  present  the  directivity  concept  in  a  manner  which  will  be 
es sen t ial  in  later  sect  ions. 

\s  a  formal  and  final  remark  we  note  that  once  H  had  been  chosen  to 
bi'  a  column  matrix,  it  was  quite  natural  to  choose  R  to  be  a  row  matrix, 
because  it  is  desirable  to  handle  only  vectors  of  f  lie  same  species,  however, 
we  can  choose  the  filter  vector  also  to  be  a  column  vector.  With  the  con¬ 
vention  that  the  frequency  responses  appearing  in  the  filter  column  vector 
are  taken  to  be  the  complex  conjugates  of  the  actual  frequency  responses, 
the  transpose  conjugate  of  the  filter  column  vector  will  be  identical  to  t  lie 
low  vector  with  which  we  have  thus  far  dealt.  Henceforth  we  shall  deal 
oniv  with  column  vectors. 
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Chapter  3 

ADAPTIVE  SIGNAL 
PROCESSING  WITH  A 
KNOWN  WAVEFRONT 

3.1 

W'  first  have  to  comment  about  stationaiity  and  adaptivity.  Consider  a 
remote  poifit  like  source  transmitting  a  random  signal  S(f).  This  signal 
reaches  a  set  ol  A  sensors  alter  passage  through  a  propagation  medium. 
I  he  geometry  of  t  lie  array  is  known  and  can  be  referred  to  in  terms  of 
coordinates.  Usually  the  spacing  of  the  sensors  is  about  half  a  wavelength 
ol  the  highest  expected  frequency.  The  outputs  of  these  sensors  are  N 
voltages  exclusively  caused  by  S{  ( )  as  modified  by  the  medium.  So  far  no 
otlu  r  source  or  noise  interference  is  presumed  to  be  present. 

I  In1  relation  between  the  N  sensor  outputs  can  then  be  described  in 
terms  ol  a  correlation  matrix  in  the  spatio-temporal  domain,  or  of  a  cross- 
spedral  density  matrix  (CSDM)  in  the  spatio-spectral  domain.  Both  de- 
-criptions  are  equivalent:  the  second  will  lie  preferred  here  for  its  simpler 
mathematical  symbolism.  Although  the  CSDM  is  used,  we  are  primarily 
interested  in  large  band  widt  hs  and  not  pure  sine  aves.  From  a  strict  point 
ol  view,  correlation  implies  stationaiity  and  neither  the  signal  S(f)  nor  the 
medium  usually  meet  this  criterion.  Nevertheless  the  limited  time  <iuasi- 
st  ii  I  mini  nt  i/  assumption  is  the  cornerstone  of  adaptive  processing  in  both 
patio  temporal  and  spatio-spectral  domains. 
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More  precisely  we  si i n  11  assume  that  the  lluot nations  of  the  medium  are 
slow  ('noiiu.li  lo  I  >e  considered  as  quasi-st  at  ionary  over  a  linn'  /.  which  [or 
I  Ik-  ocean  is  estimated  to  he  of  the  order  of  one  second,  Concerning  the 
signal  S  ( /  i .  we  also  assume  tin'  same  t  ime  T  of  quasi-st  at  ionarity.  Furl  lter- 
niort'  wr  shall  In'  concerned  with  signals  of  more  than  a  thousand  Hertz 
bandwidth  11  .  With  such  a  large  I  ime-bandwidt  h  product  (II  /  •  1000). 

the  problem  of  estimating  the  second  order  statistics  should  not  lx  too  crih- 
,\ tivwav  such  estimation  methods  have  been  t lie  focus  of  much  scient ilic 
literature  jlj  j  I  if'  and  improvements  can  still  be  reasonably  expected. 

The  present  paper  assumes  the  validity  of  these  estimation  methods, 
f'nder  these  conditions  tin'  elements  of  a  ('Sl)M  can  be  estimated  with 
little  significant  error.  It  is  then  pointless,  and  frequently  unwise,  to  im¬ 
pose'  an  artificial  relationship,  such  as  the'  Toeplitz  assumption,  upon  the 
experimental  values.  This  position  being  adopted,  two  remarks  must  be 
made. 


1.  All  the  estimation  methods  commonly  in  use  today  implicitly  rely 
on  f  lie  i  .risli  oci  oj  o  slolioiwrii  process  which  is  (  shmntcil  during  a 
limited  time  interval.  Actually  the  real  problem  is  Imw  best  to  lit 
this  stationary  process  to  a  non-sl  at  ionary  data  sequence.  Although 
this  distinction  may  seem  slight,  it  could  be  the  basis  for  refining  the 
est  imal  ion  met  hods. 

2.  Kverv  estimation  method  relies  on  assumptions  about  the  process 
influencing  the  data  outside  the  available  sequence.  It  seems  reason¬ 
able  that  a  "good"  met  hod  should  be  one  wit  h  few  assumpt  ions.  This 
concept  will  be  used  as  a  guideline  in  this  paper. 


In  short .  we  will  ohly  consider  spectral  and  cross-sped  ral  densities  which 
are  valid  for  a  limited  time  and  experimentally  renewed  every  'I  seconds  or 
last  er.  such  as  the  I )  KITIC  (  Delay  l.ine  Time  (  ’(impression  )  pn  iced  lire1 1 1  1 1. 


1  I  Ins  almost  fopjnll  »-n  )>r<  >t  «'l  hi  •*.  w  1 1  i  <  1 1  was  nj  mjnall  v  1 1  n  j  >1*  •  i  n  *  •  1 1 1  *  *  I  usini;  »l*‘lav  lines, 
.an  !>*•  H  nn^posril  lo<mnput»i  a  1*^*  »i  it  Inns.  M  pn>vi<!'-s  tin-  iiiaxiniiiiii  ainmuif  n|  ovrilap 
}  >i  ii  }4,  ‘fait  (maximum  a«laplivifv)  rompaftU'-  with  I  li»  sampling  jatr,  al  tin  <<»st  nj  s.nm 
,  Miuinilei  spf.rl  mult  in!i''l  ! » v  a  la<  Pu  « *1  *t .  i  i  n  ♦  *  1  fp»m  tin  film  l»aml  wi«|tli  pi'»«lu<t. 
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\\ «'  must  now  he  quite  clear  about  the  source  vector  and  wavefront.  Within 
the  framework  of  the  previous  quasi-slationarity  concept,  let  us  recall  the 
single  source  transmitting  S(t).  The  transfer  function  between  this  source 
and  the  output  of  the  n,h  sensor  can  be  represented  by  a.  linear  filter 
with  a  complex  frequency  response  <1 >„(/)  which  is  a  component  of  the 
A  -dimensional  complex  vector  <t >(/)• 

II  (/,(/)  is  the  true  spectral  density  at  the  source,  the  cross-spectral 
density  between  the  output  of  sensors  n  and  m  is  given  by  a  well  known 
expression  already  presented  in 

To  disencumber  the  notation,  the  frequency  will  henceforth  be  treated  as 
an  implicit  variable  except  when  the  contrary  is  specifically  stated.  With 
t  Ids  change  in  notation, 

<?.»»  =  (3.1 ) 

i>  the  n1’1  row,  m'1'  column  element  of  the  C'SDM  linked  to  the  source.  By 
convention  <1>  is  given  the  status  of  a  column  matrix,  then  the  C'SDM  is 

Q  =  (3.2) 

We  then  deline  the  “source  vector”  as 


S  =  d?Q, 


which  allow 


,s  the  C'SDM  to  be  written  as  a  dyadic  product  of  the  vector  S, 


Q  =  SS 1 .  (3.4) 

We  also  find  it  convenient  to  introduce  the  concept  of  “wavefront”  as 
nnrnuilizi  <1  vector  (and  also  column  matrix) 


where  t lie  scalar 


<l>  I  ( <1>  ^  <l> ) 


TO 


■>■> 


ADAPTIVE  PROCESSING 


HIM  I  ollV  iously  F^F  I.  Vector  F  is  not  the  physical  wavefront  but  only 
its  transfer  through  the  medium  and  array  sensors.  For  the  sake  ol  lexi¬ 
cal  simplicity  let  us  accept  (his  slight  discrepancy  with  orthodox  physics. 
The  wavefront  is  then  associated  with  some  particular  source  position  and 
associated  with  it  is  a  fixed  set  of  filters. 

The  spectral  density  dt  ol  the  source  is  unknown  since  it  is  the  object  of 
identification.  The  vector  <f>  (ol  filters  <t>„)  is  also  unknown  since  we  dismiss 
any  knowledge  of  t he  medium  in  this  paper.  What  can  he  measured  at  the 
output  of  tin'  sensors  are  the  components  of  the  source  vector 

S  d)  <l> 

df  |$|  F 

with  their  relative  amplitude  and  phases.  Bnt  with  the  absence  of  any 

1 

knowledge  about  the  source,  it  will  never  be  possible  to  isolate  dt~  from  the 
unknown  !<1>|.  So  introduction  of  the  wavefront  vector  is  justified,  and  the 
array  is  only  able  to  yield  an  apparent  spectral  density,  d.  given  by 

r/'  dj  | <l> |  .  (•?.(>) 

This  relation  relied s  the  ability  of  interference,  for  example  between  co¬ 
herent  mult ipal lis.  to  alter  the  spectral  density  of  the  source  as  it  is  seen 
through  the  array.  Finally,  the  expression  for  the  source  vector  is 

S  </*F  (T7) 


with  F'F  I  and  S*S  d. 

In  this  section,  the  wavefront  vector  is  supposed  to  be  known  a  priori, 
even  if  it  does  not  correspond  to  a  plane  or  spherical  wave.  The  goal 
of  ident ideation  is  the  analysis  of  the  only  accessible  spectral  density,  the 
apparent  one,  d. 


3.3 


At  this  point  we  combine  the  sensor  outputs  to  obtain  an  estimate  of  the 
a  p pa  l  ent  spec)  ral  density  using  <  <  uicepl  s  u lent  ioned  pian  n  >u sly  in  cliapl er  2. 
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This  is  relatively  easy  since  the  wavefront  vector  is  known  by  it.s  N  complex 
components  /„.  We  then  build  a  set  of  filters  with  the  complex  conjugate 
responses  f,] .  Adding  the  output  of  these  A’  parallel  filters  produces  the 
■global  filter 

V/„/i;  =  FtF=l  (3.8) 

M 

The  frequency  response  of  this  all-pass  filter  yields  the  apparent  spectral 
density  of  the  signal  for  the  purposes  of  analysis  and  identification.  This 
will  be  referred  to  as  the  “main  output." 

Incidentally  the  filters  /„  may  not  be  causal.  When  this  occurs  they 
must  be  built  with  the  same  additional  time  delay,  in  order  to  translate  all 
tin-  impulse  responses  equally  within  the  domain  t  >  0.  Also,  when  a.  set 
of  filters  is  claimed  to  lie  represented  by  a  column  matrix  F,  it  is,  in  fact, 

i  lie  row  matrix  with  complex  conjugate  frequency  components  of  F. 
With  this  convention,  it  is  not  necessary  (as  in  chapter  2)  to  use  different, 
symbols  for  a  wavefront  and  a  set  of  filters  matched  to  that  wavefront  (see 
sect  ion2.b ). 

Dropping  the  fiction  of  a  single  source,  the  presence  of  some  unwanted 
noise  in  the  main  output  will  now  be  considered.  Contrary  to  chapter  2, 
the  spatial  structure  of  this  noise  is  unknown.  Indeed,  it  docs  not  need  to  be 
known.  The  only  requirement  a  source  must  meet  to  be  considered  noise  is  a 
lack  of  correlation  with  S(t).  In  general,  the  source  may  be  either  spatially 
concentrated,  such  as  a  jammer  with  an  allowable  wavefront  vector  which  is 
diflerent  from  F.  or  diffuse  such  as  background  shipping  without  a  distinct 
wavefront.  T  he  structure  of  this  global  noise  will  be  adaptively  learned  by 
measures  on  the  sensor  outputs  and  estimation  of  the  CSDM.  While  the 
set  of  filters  F  is  not  adaptive,  everything  else  from  now  on  will  be,  with 
t  lie  t  ime  constant  T . 


ft 

V  * 


3.4 


I'n fortunately,  as  in  chapter  2.  we  can  never  know  whether  the  estimated 
('SI  )\I  corresponds  to  signal  plus  noise  or  only  to  noise.  Nevertheless  it  is 
possible  to  get  an  indication  of  the  noise  by  using  knowledge  of  F.  Since 
the  set  o|  fillers  with  components  /„  represented  by  a  vector  L  orthogonal 


i 


F  eliminates  t  lie  signal. 


hi  equivalent ly 


LjF  0 


V  /  /' 
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t  lie  out  put  of  the  combination  of  sensors  represented  by  L  will  be  a  •■signal- 
free  reference"  ISKU). 

Within  the  A -dimensional  space  oi  the  complex  column  vectors  which 
span  the  ('Sl)M.  the  entire  (A  -  l  (-dimensional  subspace  orthogonal  to  F 
is  available  to  vectors  of  the  L  type.  The  lull  exploitation  of  signal-free 
informal iotris  then  dependent  on  this  subspace  which  may  be  spanned  by 
a  set  ol  linearly  independent  vectors  L,.  (1  ■=.  j  •  A  I)  orthogonal  to 

F.  Since  we  may  arbitrarily  choose  tin'  filters  subject  to  the  constraint 


lJf  -  0 

we  sided  a  set  which  satisfies  the  condition  of  inter-orthogonality 


(TI  I  ) 


L  L,  -  0.  j  f-  k. 


(•'>•12) 


Additionally  we  can  further  restrict  the  filter  vectors  by  requiring  imnmv- 
rilii(ioi)  between  the  output  of  combination  L,  and  that  of  combination 

u, 


l)cl,  - ■■  0.  j  r  k. 


(3.13) 


This  point  is  discussed  in  more  detail  in  appendix  I). 

I  lie  vectors  L  must  be  derived  Irom  ( 3.  I  1  )  (3.12)  and  ( 3. 1 3 ).  One  might 
argue  that  C.  the  ('SUM  of  the  noise  alone,  is  not  experimentally  estimable. 
However  because  of  (3.1  1  ),  (3.13)  can  be  rewritten  as 


l]  C  •  (/FF*1  Lc  0.  j  /  k\ 


(3.1  I) 


where  the  matrix  in  brackets  is  a  sifz.ii ;i I  plus  noise  (\S|).M.  As  a  result  we 
are  allowed  to  introduce  the  ('SI).M  estimated  on  the  sensors  into  (3.13). 
m  f/i//f//f  .s.«s  <>/. oi/nii/  prisi/icf  or  <i  //si  i/ci . 

- <  eiil  r.u  v  In  <ni  i >•  n t  assumptions.  f  here  o-rlaiub  is  never  a  noise  alone  siluation  if 
do-  •  iii  1 1  >ii  I  of  I  In-  F  <  oinl>ination.  This  is  the  usual  pi.  liranunt  in  I<uir  ran  Re  passive 
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Let  us  now  see  under  what  conditions  (3.11),  (3.12)  and  (3.13)  can  yield 
;i  unique  solution  for  the  filter  coefficients  L j.  None  of  these  equations  can 
determine  the  arbitrary  complex  scale  factor  of  Lj.  So  the  actual  number 
of  unknown  scalars  to  derive  is  only  N  —  1.  Since  there  are  N  —  l  vectors 
L,  subject  to  this  system  of  equations,  there  are  (Ar  —  1  )2  complex  scalars 
to  lie  derived. 

Now  (3. 1 1  )  corresponds  to  N  —  1  scalar  equations  while  (3.12)  and  (3.13) 
each  correspond  to  (N  —  l)(Ar  —  2)/2  scalar  equations  obtained  as  j  steps 
from  1  to  Ar  —  1  ami  k  steps  from  1  to  j  —  1.  Totaling  the  number  of  scalar 
equations  we  get  ( Ar  —  1)  +  2(Ar  —  1 )( A^  —  2 )/ 2  =  (N  —  l)2  which  is  the 
number  of  equations  required  to  determine  all  the  complex  scalars  uniquely. 

We  then  have  a  saturated  system  of  equations  from  which  the  vectors 
L,  can  l>e  derived.  Although  these  filter  coefficients  are  completely  inde¬ 
pendent  of  signal  presence  or  absence,  they  are  adaptive  to  the  noise.  Thus 
il  is  possible  to  transform  the  original  set  of  N  sensor  outputs  to  another 
system  for  which  the  “main  output”  lets  signal  through  as  an  all-pass  fil¬ 
ter  (unaltered  signal  with  the  spectral  density  d),  and  the  other  (N  —  1) 
outputs  are  uncorrelated  with  the  signal  and  each  other  but  possibly  cor¬ 
related  with  the  noise  in  the  main  output.  From  now  on  C  will  denote  this 
1 1  ansformed  C'SDM. 


I  he  above  possibility  is  so  important  that  it  is  worth  further  illustration 
by  a  different  approach.  Let  us  first  require  the  vectors  Lj  to  be  orthogonal 
only  to  F  but  not  to  each  other,  and  to  be  distinct  so  a  full  description 
of  the  noise  subspace  is  possible.  In  particular  we  can  easily  build  the 
non -adaptive  vectors: 


j'1'  row 


(3.15) 
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which  arc  distinct  since  they  use  different  pairs  of  sensors  and  which  satisfy 

lJf  n. 

These  L,  combinations  yield  correlated  SFRs,  however  we  wish  to  retain 
the  signal  free  information  in  a  se*  of  uncorrelated  SFRs.  To  do  tins  we 
linearly  transform  a  set  of  correlated  SI‘Rs  to  a  set  of  uucorrelaled  Sh  Its 
which  are  their  “orthogonal  images.  The  method  adaptively  estimates 
the  eigenvectors  of  the  noise  (’SI)M  suhspace  orthogonal  to  F.  Then,  ev¬ 
en-  eigenvector  is  used  to  form  a  linear  combination  of  the  original  SFR>. 
Because  of  the  particular  properties  of  eigenvectors,  the  new  Sh'lts  are  no 
longer  correlated.  This  process  is  the  physical  expression  of  the  diagonal- 
ization  of  a  matrix  which  is  discussed  in  appendix  1). 

With  respect  to  this  new  set  of  basis  vectors,  the  OSDM  lias  a  specific 
aspect  involving  diagonal  elements  which  represent  spectral  densities,  and 
oil-diagonal  elements  in  the  first  row  and  column  which  represent  cross- 
spectral  densities  with  the  main  output;  all  other  oil-diagonal  elements  are 
zero  because  the  SFRs  have  been  mutually  decorrelated. 
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with  any  linear  combination  ol  noises  such  as  a  SFR.  Any  cross-spectral 
density  between  two  SFRs  is  zero  since  the  noise  has  been  decorrelated. 

W  e  reach  the  situation  illustrated  in  figure  if.  I  in  which  a  main  output 
with  a  possible  signal  of  spectra!  density  </  and  a  noise  /i()  with  a  spectral 
density  c„„  is  combined  to  produce  the  final  output  with  the  uncorrelated 
SFR's.  li |  to  liy  having  spectral  densities  en  to  I  lie  spectral  den 

sities  c„„  can  be  adaptively  measured  on  the  corresponding  out  puts  except 
e which  may  be  "polluted  by  the  signal  presence. 

We  now  build  a  vector  V  representing  a  set  ol  filters  which  is  an  ail  pass 
filter  for  t  he  signal  and  which  yields  a  minimized  noise  sped  ral  density.  I  his 
vector  has  e„  I  for  the  liist  component  and  combinations  ol  the  main 
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Figure  3.1:  Simplified  Processor  with  (N  —  l)  Decorrelatecl  Signal-Free 
References. 


ut  put  and  SFKs  for  the  remaining  components, 


V1  =[1,1'!.  e2, . . .  i'N' 


(3.17) 


The  resulting  spectral  density  projected  on  V  is  given  by  the  quadratic  form 

V  1  CV.  The  constraint  rn  —  I  on  vector  V  can  be  expressed  as  follows.  Let 
ns  consider  an  A  -dimensional  vector  II  with  all  its  components  null  except 


the  first  one.  I!  ~  [  1 .0.0. .  .  .01.  The  constraint  can  then  be  written  as 


flty  =  ytn  =  1, 


eh  are  particular  cases  of 


V1  II  —  1 . 


\\e  can  then  vary  vector  V  so  as  to  minimize  V'CV  which  is  equivalent 
to  maximizing 


V'CV 


(3.19) 


•_'S 

.1  DAPTIYi:  PttOCEPSlSi, 

The  solution  is  derived  as 

in  section  2.(> 

.  wit  h: 

X 

replaced  In- 

V. 

replaced  by 

II. 

Y 

replaced  by- 

MV. 

Po 

replaced  by 

(m* )  '  ii. 

Equation  1.1  1  can  be  written  as 

CV  -  O 

li.  i:{.2iH 

where  0  is  an  arbitrary  sc 

alar.  What  matters  is  that  (.‘5.20)  expresses  the 

fact  flial  the  components  of  CV  are  nullified  except  for  f  lie  first  one.  Fsing 
(.'{.Hi)  and  (.{.IT)  we  now  write. 

‘Vo  t  <•/',,  <>•  J  I-*'-  '-{/Jl  , 

Two  remarks  must  he  made:  first,  the  spectral  density  c,l0  <|o('s  not  ap¬ 
pear  in  (2.21  ).  hence  it  is  also  absent  from  the  components  e,.  Thus  vector 
V  represents  signal-free  information.  Second,  every  component  of  V  ( <■  \ 
cept  Co  I  )  is  dependent  only  on  the  spectral  density  ol  the  corresponding 
SFR  and  the  cross-spectral  density  of  this  SFR  with  the  noise  on  tin-  main 
output.  This  last  remark  is  intuitively  predictable  from  the  11011-correlation 
of  the  SI’ |{s.  The  entire  process  is  <1  repetition  of  independent  pronssi*  lor 
each  SFR.  The  filter  to  be  put  on  the  SFR  is  the  complex  conjugate  of 
t  he  component  c,. 


Although  we  might  now  give  a  formal  expression  of  the  minimized  noise 
V^CV  with  the  vector  V  we  have  just  determined,  it  seems  more  fruitful 
to  look  into  the  physical  significance  of  such  a  result. 


3.6 


Let  us  consider  separatelv  I  lie  mam  out  put  and  the  first  Sk  R .  I‘  iltcring  I  In- 
latter  with  c„|  n, ,  and  adding  tin-  result  to  tin-  main  output  is  <-<piival<  ret 
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MAIN  OUTPUT  WITH  B0  (sd  :  c00) 
AND  SIGNAL  (sd:d) 


B,  (sd  :  cn)  B' 


Bi 


Bs  ■+■ 
SIGNAL 


Figure  3.2:  Processor  with  a  Single  SFR  Noise  Subtracter. 


to  filtering  the  SFR  with  Coi/cu  and  subtracting  the  result  (see  figure  3.2). 
f.d  us  now  proceed  to  design  a  filter,  .r  (also  shown  in  figure  3.2),  which 
minimizes  the  spectral  density  of  the  noise  in  the  output.,  B,  =  Do  —  B[. 
Hence  using  the  rule  given  in  appendix  A  with  the  notation  sd  =  spectral 
density  and  csd  cross-spectral  density: 


sd|/*,l  csd  I  B,B,  | 

csd  |(  B0  -  B[)(Bq  -  B[ )] 

sd  \B0\  cs  (\\B\Bo)  -  cs<l|£o0!)  +  sd[^;] 


Hut: 


So  that 


csfl  [/?( Bo]  =  c10.r 

csd  [£(,#!!  =  e0i.r* 

sdf^iJ  =  c, i  | .r | 2  =  cn.r:r' 

sd  [  B, ;  c„{)  c |  „.r  r„ ,  .r '  t  c , ,  .?•  .»•  (3. 23  ) 


\’arving  .r  to  frrt  a  minimum,  we  must  have  for  any  value  of  (he  increment 


,\i).\rrivk  ritoi  ' 


hence 

which  again  justifies 


cnr  <'oi 


ll  is  remarkable  t  li;it  the  left  side  of  (3.24)  is  precisely  the  cross  sped  rai 
density  of  Hf  with  If.  Indeed  looking  at  figure  3.2: 

It,  Ho  H\ 

csd  i  B,  If  i  esd  |/?o^il  csd  |  If  If | 
csd  |  Hu  If  i  <01 

csd|/r,/M  (csdizi,y;i)' 

(<•,,.<')' 

-  cn.r 

hence 

csd  |  ffl.f\  -  <’oi  <’i  i •<' 

which  is  precisely  the  quantity  nullified  in  (3.24)  1<>  minimize  the'  output 
noise.  If  so  happens  that  tin  sunn  filter  trliich  proriihs  tin  minimum  <>/ 
output  noisi  .  also  nullifies  tlx  eorrelat  ion  hi  tin  in  tin  output  uiut  tin  Si'll 
in  rolrt  il. 

In  very  qualitative  (even  hardly  scientific)  terms  a  physical  interpreta¬ 
tion  could  he  presented  as  lollows.  I  he  filter  <'oi/<'ii  draws  from  the  SIR 
what  could  he  considered  as  the  hest  "likeness  of  II n  so  that  the  suhtrac- 
tjon  of  this  likeness  yields  a  minimum  of  residual  noise.  As  a  result,  the 
"likeness”  of  the  SFR  has  hern  removed  from  the  residual  noise  remaining 
in  II,,  in  other  words  II,  finds  itself  decorrelated  from  the  SKK. 

'The  spectral  density  of  B,  as  derived  from  (3.23)  and  (3.21)  is 


|c„|  1 

Cl  1 

( 3.23 

r  I  I  .  We  can  see  1  ll 

at  lot  t  his  par 

rl iin/  1  iro  noiM  s 

IS  t  <j  If  I  I'd  hill  1 1 

far  from  a  genera 

1  rule  a  ltd  Ilia  c 

at  first  glance.  e\<  n  he  looked  upon  as  paradoxical. 
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3.7 


I  hr  repetition  of  (  tie  same  process  with  I  lie  other  .SFRs  is  obvious.  With 
/  stepping  from  I  to  ( A'  -  1),  we  put  a  filter,  coj/cjj,  on  the  jtU  SFR, 
IF.  ami  subtract  the  result  Z?'  from  the  main  output  to  obtain  the  set  of 
spectral  density  reductions.  |c0j|"  /('jj-  When  all  the  SFRs  are  exhausted, 
1  he  out  put  of  t  he  whole  process  is  uncorrelated  with  each  SFR  and  possesses 
the  minimized  spectral  density 


C min  “  '  00 


V  m! 

,(=i  rn 


(3.26) 


It  is  easy  to  check  that  rm;n  is  positive  and  precisely  the  value  of  V^CV 
(from  section  3. a)  when  V  is  the  solution  of  (3.22).  Since  the  signal  is  not 
correlated  with  the  SFRs.  the  unaltered  signal  spectral  density  d  may  be 
addil ivrly  superposed  on  cm;„.  While  the  sum  on  the  right  side  of  (3.26) 
only  involves  signal  free  estimates,  c0o  may  or  may  not  contain  signal.  As 
a  result  we  still  can  not  say  for  sure  whether  the  spectral  analysis  of  the 
linn)  output  yields  cm;„  4-  d  or  c„„„  alone. 

Although  we  are  in  a  common  situation  with  a  signal  and  a  minimized 
noise  under  the  constraint  of  unaltered  signal  (Wiener),  we  are  close  to 
achieving  something  better.  Indeed  only  one  scalar  relation  between  the 
elements  of  the  noise  alone  matrix  C  woidd  provide  a  way  of  computing 
noise  alone.  r,m.  as  a  timet  ion  of  t  he  ot her  elements.  We  have  exhausted  the 
lull  array  power  and  now  need  some  a  priori  information.  Incidentally,  it  is 
easy  to  check  that  the  wanted  relation  is  not  provided  by  the  computation  of 
o  and  of  the  first  component  <>t  vector  CV  in  (3.20).  Soil  any  bit  of  reliable 
information  is  available  about  the  original  noise,  it  can  be  transposed  to 
the  SFRs  and  provide  a  relation.  That  would  give  us  a  “super-resolution” 
system  where  the  second  order  statistics  of  the  unwanted  noise  could  be 
eliminated,  and  the  spectral  density  of  the  signal  could  be  derived  from  the 
spectral  analysis  of  the  final  output  with,  in  principle,  no  error. 


3.8 


If  no  assumption  is  suggested  by  the  physics  of  the  situation,  we  can  still 
tr\  to  illustrate  how  the  “super-resolution  system,  yielding  no  error  in 


1 1>  \  r !  1 1  i  nh)>  /  '-•>/ 


>11411. il  ii  lent  i  |i<  a 1 1>  >11 .  could  lx-  e'onsiilereel  an  "as  \  1 1 1 1  >  1 1  >t  n  a  1 1  \  large-  aim. 
assumption."  Tin-  original  CSDM.  I',  t  a  k*n  directly  lieliiiul  I  In-  sensor-  in 
,i  Inilv  operational  situation  cannot  show  an  i.ni<lli/  mill  e'igenvalm.  Hut 
t  lie  smallest  initv  lie  very  small  a>  compared  to  the  largest  mu  .  In  la<  I  .. 
tost  that  t  lie  ari  a v  has  enough  sensors  to  "surmount  t  lie  complexity  o|  i ! «« 
acoustic  lield  should  appear  when  looking  at  the  set  of  eigenvalues.  I  in¬ 
set  should  exhibit  a  large  span  of  values.  Intuitively  the  smallest  should 
decrease  wlu'ii  the  number  <>l  sensors.  A.  increases.  I  lie  validity  ol  this 
"asvmptot  icallv  lai  ge  array  assumption  is  dependent  not  <ud\  on  \  bene.: 
m  uch  la  rger  than  the  n  umber  ol  sources,  but  also  on  the  u  iedi  u  u  i  cm  1 1  ph '  .  1 1  v 
and  the  struct  lire  of  the  dill  use  noise  sources  (  s  ee  chapter  I  i. 

Now  for  any  ('SI)M  with  a  non-zero  smallest  eigenvalue  we  can 

write' 

I  I  1  inm^.Y  * 

where  I\  stands  lor  the  unit  matrix  ol  order  A.  Then  I  is  a  "singular 
matrix  with  one  zero  eigenvalue’  or  more'.  A  we'll  known  plivshal  interpre1 
lation  ean  be  preipeiseel  in  which  's  ' he  spectral  density  <  >1  a  Irae-tion 
of  noise'  which  is  uncorre-lateel  betwe'e-n  sense »rs  and  has  the’  same'  speed  ral 
densitv  on  eaeli  se>nse»r.  Ibis  i  n  t  e  ■  r  p  r<'t  at  ion  e  >f  emill  is  typieai  e.|  a  sort  of 
bae  kgreninel  imise  of  separate  but  identical  pre'ampliliers.  Hut  such  an  «  le  e- 
I  l  ie  noise'  is  known  le.  lie-  far  be-hiw  the  minimum  acemslic  tewed  in  nm.  •  ! 
arrays.  Sueh  an  inlerpre't  at  ion  <>l  eleies  not  m<'an  that  it  is  t  h »  /i/ii/wcei/ 

n  a  lil  il  for  any  CSDM.  I.  but  as  the  array  grows  in  size'  and  the  span  e>|  llo 
e'igi  n value's  spre’aels.  it  be-conu's  more  ainl  1 1 ie >r<-  like-iy  that  tin*  edee-nun-mg 
,P1]M  lolllel  re-ally  represe-llt  this  kinel  e»f  ne>is«'.  It  Weillld  lie'  still  llion-  ele'eli 
|>|i'  if  t  we i  or  l  hive  eigenvalue's  were  found  with  almost  the*  same-  iiiiiiiuiuin 
value':  ami  eel  eeuirse.  a  minimum  prae  tie-ally  ne'gligible'  as  compare-el  to  tin 
largest  e’ige'iivalne’s.  The  "asympteitieallv  large-  array  assum|itiou  eeuisists 
in  stating  that  Hits  sort  of  noi.tr  cunnat  <  .ns!  uroii.t/inilh/  and  therefore'  the- 
v i e •  1  e  1  e •  e l  b\  the-  original  array  ean  be-  simply  m-ghe-t  e-el . 

In  tin-  frame  e»f  sueh  an  assumption  1  lies  <nm's.  like'  I  .  e  "singular 
(  ’S  I )  M  .  Since'  tin-  main  e  >u  t  pu  t  ami  t  h  <  ■  N  f  K  -  are  'hri  vee  I  lim  a  rl  \  I  coin  tin 
-i  lion  output*-,  tin  <  S  | )  M  C  ■  i>  a  bo  simgula  i  .  w  Im  h  iman- 

< I e ■ t  I  C  -  i  It.  1  ■  •  ’  i  i 

•  )  i  <  r  >  "del  "  d<  1 1  •  -t  e  -  t  In'  <  1  •  •  I  e-r  m  ilia  lit  <>|  I  lu  matrix  1  ■  •  a<  •  e  pt  I  .1  e  i  e- 
to  .leu-  ilm  e-,i  i.-n'i  ol  t  hi-  wr.  si  i  i  n  v; « ■  u  1 1  ■  de-tf-i  iiiiii''il  "luukgioeind 
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aeonst  ic-noise  on  the  main  output  and  the  SFRs.  It  is  easy  to  check  that 
I  .'l.'Jti)  also  reads 

<!<*<(  C)  . 


I’l  1 ('22  •  •  -  jY ' .V ' 


(3.2*) 


so  that  accepting  (3,. 27)  is  simply  accepting  that  there  is  no  noise  on  the 
linal  output  and  we  have  got  a  “super-resolution"  system.  This  can  only  be 
credible  as  a  limiting  representation  of  infinitesimal  “background”  acoustic 
noise  if  C  possesses  a  large  span  of  eigenvalues  and  if  the  minimum  value 
is  reasonable  negligible. 


I, el  ns  now  reconsider  figure  3.2  which  shows  a  “noise  subt  ractor.”  This  ter¬ 
minology  is  generally  associated  with  “noise  canceling”  techniques,  where 
the  unwanted  noise  on  the  final  output  is  really  nullified.  In  fact  noise  can¬ 
celing  is  the  same  device  as  in  figure  3.2  but  used  in  a  particularly  favorable 
situation  where  B\  and  B n  are  “fully  correlated.”  It  means  that  both  noises 
coming  ii<  m  one  source  through  linear  filtering,  B0  can  be  seen  as  a  noise 
produced  from  B\  through  a  linear  filter  <A.  We  then  have 


csdi/7,/?,,! 
sd  I  Bo) 


('11  M2 


and  t  herefore 


roor  1 1  '  kiol  (  3.20 ) 

for  any  possible  o.  This  is  the  usual  test  of  “full  correlation.”  Equation 
I  3.2b)  then  yields  sdj/7,,j  ;  0  and  the  noise  subtractor  is  as  nearly  perfect 
as  lull  correlation  is  achieved. 

It  is  time  now  to  complete  figure  3.2  with  the  processors  it  implies. 
This  is  done  in  figure  3.3.  The  two  time  delay  devices  (or  algorithms)  are 
there  to  ensure  that  the  subtraction  is  really  achieved  with  the  very  same 
sequences  of  /7n  and  /7 ,  which  participated  in  the  building  of  the  filter  with 
a  time  constant  7  .  This  device,  or  the  algorithmic  equivalent  of  figure  3.3, 
is  known  in  the  French  scient ific  literal ure  as  a  “correlo- filt  re  and  has  been 
extensively  studied  and  improved  in  the  past  fifteen  years  |Jb-17). 

\  particular  feature  of  the  correlo-filter.  in  fact  a  degree  o'  freedom  for 
l  In-  processor,  is  that  any  preliminary  linear  filtering  of  the  SKK  does  not 
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change  the  spectral  density  of  H,.  This  is  why.  in  its  present  form,  we  have 
pirderred  to  Jirsl  "self-whiten''  the  SFK.  H |.  with  an  open  or  closed  loop 
algoril  lim  and  then  proceed  wit  h  a  normalized  value  of  I  lie  spec  I  ra  I  density, 
rn  e.  which  is  not  fiiqtinuti  <l<  j><  nth  nl  in  the  uselul  liandwidt  h. 

A  self- whitening  filter  must  have  a  frequency  response  II,  such  that 
dl',!2  e  r  1 1 .  The  phase  law  of  this  filter  is  free  which  makes  it  easier  to 
Imild.  and.  as  mentioned  above,  does  not  alleel  I  In'  result.  1‘igure  is 

then  modified  as  shown  in  limine  If  I.  where  the  filtering  of  I  lie  whitened 
SI  H  ll\  can  lie  conceived  in  the  time  domain  as  its  convolution  with  the 
cross-correlat ion  function  of  //„  with  U\  and  the  essential  scale  laelor  I  e. 

If  all  the  other  Sf'lts  are  whitened  at  I  lie  same  spectral  density  c  and  it 
we  name  correlo- filter  ((•!*')  that  which  is  f  ramed  by  t  he  dot  ted  line  in  figure 
I.  the  processing  from  the  sensor  out  put  s  can  be  sketched  as  in  figure  TV 
This  figure  could  be  presented  as  an  open  loop  solut ion'  to  the  original 
problem  of  gel  t  ing  t  he  signal  sped  ml  den  si  I  y  with  a  minimum  of  mi  want  ed 
noise  and.  for  a  large  array,  wit  li  prad  ically  im  noise  at  all.  Neverl  lieless  it  is 
not  intended  here  to  propose  such  an  open  loop  algorithm  for  operational 
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amplitude  signal  superposed  <mi  /?,  which  can  l.e  considered  as  derived  from 
tin1  main  signal  throimli  a  Idler  whose  frequency  response  will  lie  denoted 
I’.v  I  hen  We  have  a  new  csd  between  the  main  output  and  the  SI  R. 

Sn  ‘  in  •  <!'■'  •  ( •{.■!()  I 

while  I  he  spectral  density  is  also  modified, 

Si  Si  !  ^  )*-')'  .  (d.d!) 

I  oyet  tier  and  (  d.d  I  )  are  the  < pi  ant  it  ies  est  imal  ed  ||\  t  lie  correlo  lilt  e| 

w  1 1  i ** 1 1  y  adds  a  lilt  er  (J  <  n  |  e ,  | .  I  he  lorn  n  r  on  I  pn  I  noise  /PJ  is  a  1 1  erec  I  |  ( . 
Ii\  so  that 


s,l  f /• 

L 

f,  n\] 

sd  j  /*,]  (  > 

<  uO  . 

il  IK 

sd 

[/>■; 

-d  |  /?, |  <  >  - 

>'i  i  (J  . 
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On  1  hf  main  output.  B  becomes  B,  with: 


sd  [/?,]  =  scl  B0)  -  cscl  [B0B;]  -  cscl  [B\B0]  f  sd  [/?;]  , 
csd[B0fi|]  ^csdffioBjQ-  =  c0lQ-, 
sd  !i?s]  —  Cqo  ~  c-oiQ  ~  <\uQTciil4?l  • 


Irnrr 


sd  Bs  ~  c oo  — 


b’o.i 


Oi 


and  1  i  nail  v 


sd  (  B,  I  —  Cqo 


Kot 


+  </((/’  +  r )  +  d2  iv-is 


(3.32) 


<’ii  +  d|0|" 

noise  itself  is  increased  by  only  a  second  order  term  when  0  is  small, 


<'oo  ~ 


koi  I 


1 1 


+  d  1 0 1 


/.  1 2  ‘ 


The  main  point  is  that  some  signal  now  comes  through  the  loop  of  the 
correlo-lilter.  This  produces  an  additive  alteration  of  the  spectral-density 
with  the  first  order  component  d(  0  f  0')/e  n.  So  the  signal  spectral 
density,  dll  -  (c  -t-  f')/cu|.  may  be  larger  or  smaller  than  the  expected 
density,  d.  in  an  unpredictable  way  since  nothing  is  generally  known  about 
the  filter  Small  or  large,  the  most  spoiling  effect.  of  signal  leakage  into  an 
SI  If  lies  in  the  fact  that  the  signal  spectral  density  becomes  distorted  by  a 
process  which  has  been  conceived  to  preserve  it  .  That  is  why  other  met  hods 
have  to  be  considered  when  the  signal  wavefront  is  unknown  (chapter  4). 


3.11 

Nevertheless  compromises  can  be  considered  in  order  to  make  the  previous 
processing  more  robust.  Leakage  is  a  consequence  of  uncertainties  about 
the  signal  wavefront  which  can  be  expressed  in  terms  of  subspaces  and 
ST  Its.  If  the  signal  wavefront  cannot  be  exactly  assigned  the  shape  of 
vector  F.  at  least  it  will  be  generally  possible  to  accept  it  as  being  located 
in  a  limited  subspace  “around"  F.  This  signal  subspace  is  the  genera lizat ion 
of  the  notion  of  beamwidth  inasmuch  as  the  latter  concerns  the  concept  of 
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aperture  limitation  and  reflects  only  the  uncertainty  about  the  direction 
of  an  incoming  plane- wave  signal.  Anyway,  the  signal  siibspace  lias  to 

lie  determined  by  a  set  ol  vectors  {F[,Fn . F(,},  />  ■■  A.  more  or  h’ss 

neighborin';  F. 

The  signal-free  complementary  subspace  is  orthogonal  to  all  these  vec¬ 
tors.  and  is  consequently  an  (A  />)-dimensionnl  subspace.  With  only 
I  A  /< )  sensor  combinations  being  considered  as  SFRs  the  process  01  re¬ 
moving  I  lie  unwanted  noise  Ironi  all  t  he  />  “main-out  put  s  '  is  likely  to  be  less 
ellicient  than  Indore.  Indeed  the  reduction  ol  the  SFR  number  (for  a  given 
A  )  lessens  the  possibility  of  exploiting  an  ideal  “quasi-singular"  CSDM 
of  unwanted  noise  and  reaching  the  threshold  where  only  a  questionable 
acoustic  “background"  noise  could  not  be  eliminated.  Such  is  I  In'  cost  to 
be  paid  for  the  uncertainty,  plus  the  painful  repetition  ol  the  noise  subtrac¬ 
tion  process  from  the  />  main  outputs  which  requires  a  total  ol  p{  A  />) 
correlo-lilters  instead  ol  (A  I). 

In  short  (he  larger  the  uncertainty,  the  larger  the  signal  subspace  and 
the  smaller  the  signal  free  subspace  whose  size  is  a  measure  of  the  process 
efficiency.  Total  ignorance  of  vector  F  would  lead  to  a  generalization  ol 
beam  forming  in  an  A' ’-dimensional  complex  space:  a  mosaic  of  non  overlap¬ 
ping  signal  subspaces  each  one  being  processed  after  its  own  complemen¬ 
tary  one.  |,et  us  remember  we  are  considering  large  arrays  ol  (it)  sensors 
or  more.  The  computational  effort  seems  beyond  imagination,  a  true  over¬ 
dose  of  processing.  The  number  ol  sources  possessing  a  wavelront  has  to 
lie  much  smaller  than  A  in  an  array  really  able  to  cope  with  the  source 
and  medium  complexity.  All  of  these  problems  might  help  one  to  favorably 
consider  chapter  I  in  spite  of  its  own  perplexing  difficult ies. 


3.12 


As  a  final  comment  about  adaptive  processing  it  must  be  mentioned  that 
I  lie  t  eel 1 1 1 iq ties  of  adapt  i ve  s<'l(  noise  reduct  ion  I  or  cam  eling  )  are  only  a  pa r 
tieular  case  within  the  framework  <>l  the  presentation  given  in  this  section. 
As  another  example  consider  the  vibration  <>l  a  ship  structure  in  which  some 
of  the  induced  aeoiist  ic  energy  lies  within  the  signal  I  in  ml  w  id  t  h  and  reaches 
1  he  acoustic  sensors  of  a  huff-mounted  array.  Obviously  the  wise  tiling  to 
do  lirst.  is  to  merlin  n  nail  v  damp  the  source  o|  noise.  However  when  I  lie 
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best  possible  result  has  been  so  obtained  residual  self  noise  may  persist. 

One  possible  wav  to  try  and  clear  the  array  of  this  noise  is  to  use  ad¬ 
ditional  sensors  (accelerometers)  located  on  or  close  to  the  vibrating  struc¬ 
tures.  Such  sensors  have  two  very  interesting  properties  for  our  purposes. 
Thev  are  nearly  perfect  SFRs  because  they  are  mounted  inside  the  hull 
which  makes  them  practically  impervious  to  interference  from  acoustic  sig¬ 
nals  produced  outside  the  hull.  Also,  one  can  expect  that  each  particular 
output  is  highly  correlated  with  that  part  of  the  self  noise  generated  by  the 
1  est  eil  st  met  ure. 

Therefore  a  procedure  of  self  noise  reduction  can  be  considered.  It 
appears  just  as  an  increase  in  the  ‘'size"  of  the  array,  each  additional  sensor 
delivering  a  SFR.  Self  noise  SFRs  may  or  may  not  be  correlated  either 
between  each  other  or  with  the  SFRs  already  extracted  from  the  array.  If 
no  reliable  information  is  at  hand  the  diagonalizat ion  process  has  to  be 
applied  to  all  the  SFRs  for  them  to  become  uncorrelated.  They  then  feed 
the  subsequent  noise  subtracting  process  ol  sections  3.(3  and  3.7. 


3.13 


,\s  conclusion  to  clinplei  the  ma  jor  points  of  this  type  of  processing  are 
listed. 


I .  I  b  opping  unreliable  assumpt ions  about  unwanted  noise  leads  to  ada.p- 
tive  processing  with  the  necessary  complexity  required  to  track  the 
variable  unwanted  noise  spatial  structure.  There  is  no  practical  dif¬ 
ference  between  unwanted  noise  sources  gifted  with  a  wavefront,  and 
diffuse  noise.  So  the  ambiguity  linked  to  the  directivity  concept,  and 
mentioned  in  chapter  2,  disappears. 

The  processor  yielding  an  unaltered  signal  in  a  minimum  of  unwanted 
noise  can  be  analyzed  by  starting  with  t  lie  building  ol  a  specific  sensor 
combination  to  carry  I  lie  signal  and  (A’  1  )  signal-free  combinations 

(SFRs).  Then  the  SFRs.  generally  correlated,  are  combined  in  such 
a  way  as  to  decorrelate  them.  Further  processing  is  a  repetition  ol 
"noise  subtraction"  I  noise  itself  and  its  spectral  density)  between  the 
main  output  and  the  decorrelated  Sf  Rs. 


ADAPTIVE  PROCESS IX (! 


The  possibility  of 'getting  a  “perfect"  result  (no  additional  noise  mixed 
with  the  signal  I  is  linked  to  the  apt  it  tide  ol  the  array  to  surmount  the 
com jtlexil  v  ol  t  lit'  source-medium  system.  This  aptitude  is  obviously 
bounded  by  the  size  of  the  array  (its  own  “complexity")-  It  <an 
be  tested  with  the  set  of  eigenvalues  of  the  CSDM.  I  he  smallest 
eigenvalue  has  to  be  negligible  as  compared  to  the  largest  ones. 

'['In'  signal  source- vector  is  known;  but  therein  lies  the  performance 
limitation.  Any  slight  swerving  of  the  source  vector  has  to  be.  wln-si 
unpredirt  able,  counteracted  at  the  cost  ol  a  loss  ol  unwanted  noise 
reduction  possibilities.  Being  more  matched  to  the  unwanted  noise 
than  classical  beamforming,  the  processor  is  sensitive  to  any  strict 
assumption  about  the  signal  source  vector. 

It  is  a  vanity  to  expect,  in  real  conditions,  a  simple  way  o|  coping 
simultaneously  with  “imperfections"  ol  a  given  source  vector,  and  an 
unknown  multiplicity  of  sources.  The  transposition  of  beamforming 
in  a  A'-dimensional  space  is  fantastically  complex  for  large  arrays, 
and  useless  since  the  number  ol  sources  is  much  less  than  the  number 
of  sensors,  a  definition  in  some  way.  ol  an  “efficient"  array. 

The  processor  composed  of  correlo-lilters  described  in  figures  .T2.  TT 
.{.I  and  To.  is  one  possible  type  of  minimum  variance  distort  ionless 
tiller.  For  the  sake  of  physical  significance',  it  has  been  presented  ami 
discussed  here  in  hardware'  form.  Implementation  in  software  is  also 
possible  if  enough  computing  capacity  is  available. 


Chapter  4 

IDENTIFICATION  OF 
SIGNALS  WITH  UNKNOWN 
WAVEFRONTS 


\\ V  now  conic  to  tie  main  goal  of  this  paper.  How  can  we  identify  the 
spectral  density  of  a  source  with  a  totally  unknown  wavefront  which  must 
also  he  determined?  In  other  words  we  will  be  trying  to  find  the  minimum 
information  needed  to  restore  the  spectral  density  after  dropping  all  a  priori 
knowledge  concerning  both  the  signal  and  noise. 

In  previous  sections  the  source  had  been  characterized  by  a  known  vec¬ 
tor  involving  the  notions  of  spectral  density  and  wavefront.  Here  such  a 
specific  characterization  will  not  be  assumed.  As  a  consequence  there  is  no 
longer  any  way  to  distinguish  between  several  sources.  This  puts  us  in  the 
position  of  being  able  to  identify  either  none,  or  all  sources  simultaneously. 
Wilder  these  conditions  the  concept  of  what  may  be  a  signal  has  to  be  rede- 
fluid  to  include  every  source  gifted  with  an  a  priori  unknown  wavefront.  It 
then  becomes  difficult  to  separate  those  wavefronts  in  which  we  have  some 
practical  interest,  from  those  we  wish  to  ignore.  Previously  the  existence 
of  a  wavefront  which  characterizes  “interesting  sources”  has  been  implicit 
in  all  spatial  processing.  But  since  it  is  now  the  only  assumption  retained 
about  a  signal,  it  is  worth  mentioning  that  even  this  constitutes  a  priori 
knoirli  dip .  If  the  existence  of  wavefronts  comes  into  question  the  basis  of 
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all  spat i a  1  processing  is  remove*!,  including  this  whole  paper. 

Obviously  some  experiments  are  necessary  to  determine  how  much  we 
can  rely  on  the  wavefront  notion.  The  characteristic  feature  of  a  source 

gifted  with  a  “perfect"  wavefront  is  a  rank-one  CSDM,  SS^  regardless  of 
the  the  propagation  complexity  between  the  source  and  the  array.  The 
dyadic  properly  of  this  matrix  is  the  result  of  so  called  “perfect  spatial 
coherence"  which  is  characteristic  of  a  single  “perfect  source.” 

I  he  question  which  needs  to  be  tested  in  operational  conditions  is 
whether  or  not  a  remote  source,  which  is  powerful  enough  to  produce  a 
definitely  dominant  contribution  (  f 20  dB  or  more)  on  the  output  of  erery 
st  n. tor,  really  yields  a  rank-one  experimental  CSDM.  From  such  experi¬ 
ments  the  maximum  array  size  for  which  the  dyadic  source  approximation 
will  remain  acceptable  can  be  determined  to  some  degree  of  precision.  We 
expect  this  array  size  limitation  to  be  much  larger  than  the  sort  of  apertures 
on  which  every  wavefront  has  to  be  a  plane  or  spherical  wave.  Hopefully  we 
are  now  in  the  realm  of  “large  arrays”  which  has  already  been  mentioned 
in  chapter  •'!  and  will  be  more  precisely  presented  in  section  4..'!. 


4.2 

It  is  now  clear  precisely  what  properties  a  perfect  source  signal  has.  Noise 
is  everything  else,  a  mixture  of  every  kind  of  diffuse  source  not  possessing 
the  wavefront  notion,  with  a  CSDM  generally  exhibiting  a  rank  which  ap¬ 
proaches  its  order.  We  thus  have  to  disentangle  an  acoustic  field  made  of  an 
unknown  number  of  “perfect  sources"  from  a  “noise"  with  no  characterized 
wavefront . 

We  have  already  noticed  in  chapter  3  that  the  apparent  spectral  density, 
<1.  and  normalized  wavefront  vector,  F,  of  a  “perfect  source”  are  the  most  we 
can  expect  to  determine  using  an  array.  Since  (2.7)  is  S  -  cDF,  vectors  S 

and  F  may  both  be  represented  as  enlntiiii  matrices  such  that  F^F  I.  We 
now  assume  that  we  have  a  liinilttl  ninnlxr  of  inieorit  latt  </  perftcl  soureis 
which  produce  a  field  received  by  a  large  array.  By  limited  we  mean  the 
number  of  sources,  is  less  than  the  number  of  sensors  in  the  array. 
Iy  V.  which  clarifies  our  notion  of  a  large  array. 

1  he  CSDM  of  I  lie  perfect  sources  we  shall  be  working  with  can  then  be 
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represented  as  a  sum  of  dvads, 

C.  =  EMJ,  (4.i) 

P-1 

or  equivalently, 

a,  =  f;rfPF,F).  H.2) 

p=i 

I'lte  rank  of  t his  matrix  is  P,  since  every  vector  orthogonal  to  all  the  Fp  is 
transformed  into  a  mill  vector.  This  statement  assumes  the  set  of  vectors 
Fp  have  no  linear  dependence  which  would  reduce  the  rank  of  C,.  This 
depends  somewhat  on  propagation  conditions  and  can  occur  if  two  farfield 
sources  and  the  array  should  happen  to  lie  upon  a  straight  line.  In  practice 
this  is  not  a  stringent  assumption.  But  before  exploring  matrix  C,  we  have 
to  extract  it  from  the  raw  data  matrix,  the  data  C'SDM  C,.. 


4.3 

I  he  data  PSDM  includes  the  source  dyads  plus  the  C'SDM  of  diffuse  noise 
supposed  1  a  hr  uncorrelatcd  with  any  sources.  The  net  result  is  that  it  can¬ 
not  have  a  rank  lower  than  its  order  N.  As  already  mentioned  in  chapter  3 
we  can  never  find  an  experimental  CSDM  with  an  exactly  null  eigenvalue. 
I’lie  only  serious  accident  which  can  happen  to  a  data  C'SDM  would  be  the 
occurrence  of  a  negative  eigenvalue.  This  would  be  unacceptable  because  it 
means  t  he  correaponding  eigenvector  is  defining  a  set  of  filters  which  passes 
a  mgative  spectral  density  on  its  output.  Fortunately  this  accident  cannot 
happen  to  the  maximum  likelihood  estimate  of  a  data  C'SDM  because  such 
an  estimate  is  the  sum  of  many  data  dyads.  We  then  have  an  estimated 
noise  CSDM  with  the  right  features  required  by  physics:  Hermitian  sym¬ 
metry  and  real  positive  eigenvalues. 

It  is  time  now  to  come  back  to  the  notion  of  “large"  arrays  which  involves 
I  In'  concept  of  an  array  possessing  enough  sensors  to  surmount  the  complex¬ 
ity  of  a  system  composed  of  perfect  sources  plus  some  diffuse  noise  structure . 
We  must  emphasize  that  in  a  stable  situation,  the  smallest  eigenvalue  of 
a  data  CSDM  generally  decreases  whin  the  number  of  sensors  increases. 
bet  us  try  to  illustrate  it.  with  a  minimum  of  mathematical  symbolism. 
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Recalling  that  eigenvalues  represent  spectral  density,  it  is  well  known  that 
the  filter  combinat ion  derived  from  a  normalized  vector  which  di  lines  tin 
minimum  spirt  nil  density  is  precisely  the  comhinat  ion  derived  from  the 
normalized  eigenvector  associated  with  the  minimum  eigenvalue. 

Indeed  il  we  look  for  the  comhinat  ion  derived  from  a  normalized  vector 
v  such  that  v^v  l.  the  minimum  spectral  density  of  the  output  vd'.v 
is  given  hy  a  solution  to  the  type  of  variations  problem  previously  seen  in 
both  chapters  2  and  •!.  The  solution  is  now  well  known: 

Crv  Av 

vv..  .  a  11:11 

with  A  a  scalar.  Necessarily  v  is  an  eigenvector  and  A  is  the  corresponding 
eigenvalue.  Since  we  want  a  minimum  out  put  from  the  filler  combination. 
A  must  be  the  minimum  eigenvalue  and  v  the  corresponding  eigenvector. 

In  short  the  minimum  ei/jen  value  is  the  minimum  output  winch  cun  h< 
derived  from  u  normalized  combination  of  sensors.  Such  a  minimization 
applies  to  the  spectral  density  at  <  very  frequency.  Let  us  assume  we  have 
built  a  normalized  combination  for  an  A  sensor  array,  with  Hu  its  output 
and  furl  her  assume  we  can  obtain  an  addit  ion  a  I  (A-  i  I  )fl'  sensor  with  out  put 
/?,.  There  is  no  reason  for  B\  to  be  uncorrelated  with  Hu-  Therefore  using 
a  correlo  filter  such  as  the  one  described  in  figure  b  i  of  chapter  -'I,  with  Hu 
and  H\  at  the  same  places,  we  know  it  is  possible  to  reduce  the  spectral 
density  of  Hu  to  a  smaller  one,  H„.  Although  < lie  notions  of  known  signal 
wavefront  and  SFR  have  vanished,  the  C'F  is  not  prevented  from  playing 
its  basic  part.  We  also  remark  that  the  corrclo-filter  effect  does  not  depend 
at  all  on  the  level  of  H\.  but  only  on  the  amount  of  correlation  between  Hu 
and  /i| .  I  he  level  can  be  adjusted  to  a  value  corresponding  to  a  normalized 
(  A  )  1  (-sensor  combination  because  an  arbitrary  filter  is  allowed  at  the  SFH 
input  of  the  correlo-filter. 

d  ims  we  have  built  a  normalized  combination  of  (  N  f  1)  sensors  with  an 
output  spectra!  le  n  I  smaller  than  the  small/ st  ohiainahh  with  A  sensors. 
Audit  is  not  even  t  he  smallest  possible  one  wit  li  | A  i  I  )  sensors,  since  t  his 
is  given  by  t  he  smallest  eigen  value  of  I  he  (  A  •  I  )  sensor  ( 'SI  )\l .  In  short . 
whin  tin  numhir  of  si  nsovs  men  lists  tin  small/  st  inpnralui  dt  c  r<  asi  s,  or 
at  the  worst  remains  stable.  Adding  more  sensors  increases  I  lie  possibility 
of  having  a  large  span  of  eigenvalues  with  the  smallest  one  becoming  much 
smaller  than  the  largest  one.  \s  a  limit  it  leads  cither  to  a  singular  <  SI  >  M 
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with  a  vanishing  smallest  eigenvalue,  or  to  a  stabilization  of  the  smallest 
eigenvalues  which  are  repeated  for  each  additional  sensor.  In  both  cases  it 
is  a  test  t  hat  the  array  has  become  large  enough  to  overcome  the  complexity 
of  the  system  of  sources  and  the  noise  structure. 


4.4 


Having  clarified  what  we  mean  by  a  large  array,  we  shall  now  assume  such 
arnujs  for  the  remainder  of  this  paper.  Since  the  array  must  separate 
sources  from  the  noise  and  also  describe  the  source  wavefronts,  it  is  un¬ 
likely  that  a  simple  linear  array,  even  a  long  one,  could  cope  with  a  fine 
description  of  intricate  wavefronts.  Volumetric  arrays  are  more  desirable  for 
our  purposes  since  they  can  be  arranged  to  break  any  symmetry  which  may 
promote  linear  dependence  between  the  wavefronts  of  two  distinct  sources. 

We  now  have  to  estimate  the  number  of  sources,  and  in  the  absence 
of  any  a  priori  information  it  is  clear  that  we  will  have  to  gamble.  If  a 
set  of  several  eigenvalues  with  the  same  minimum  value  is  supporting  the 
acoustical  background  noise  (this  is  improbable)  the  difference  between  the 
data  CSDM  and  cmil,I.v. 

Cr  <mi»IjV  =  C,  (4.4) 

is  a  singular  matrix  of  rank  “A;  minus  the  number  of  equal  eigenvalues  c,,,;,,.” 
Mere  <min  is  this  minimum  eigenvalue  and  Ij\-  the  unit  matrix  of  order  N. 
It  is  then  possible  to  consider  C,  as  the  perfect  sources  matrix  of  (4.1 )  and 
(  1.2).  Unfortunately  there  is  little  chance  the  noise  structure  can  be  simply 
described  by  I,v-  Several  more  elaborate  approaches  have  been  suggested 
1 1  111  !‘23|  which  are  all  smart,  but  all  based  on  some  arbitrary  option. 

The  principle  behind  the  above  approach  is  a  parameterized  noise  CSDM 
whose  general  form  is  known  but  where  several  parameters  are  left,  free,  par¬ 
ticularly  a  scale  factor.  Such  a  parameterized  noise  CSDM  may  be  written 
as 

M  —  rrG((/.  />)  (  4  . -r* ) 

wlu  rc  er  is  a  scale  factor.  G  the  general  form  of  a  normalized  CSDM,  with 
a  and  I)  as  free  specific  parameters.  Then  in  the  difference  matrix 

C,  Cr  M.  (I.ti) 
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the  values  of  the  parameters  rr,  a  and  />,  are  varied  to  obtain  the  lowest 
possible  rank  of  C,.  In  this  particular  case  the  rank  is  A  •’>  but  il  the 
general  noise  matrix  G,  is  really  close  to  the  physics  of  t  he  situation  il  might 
well  happen  that  the  rank  of  C,  will  be  still  lower.  We  will  not  comment 
about  the  computing  dif  lien  It  ies  in  solving  (  l.(i)  since  il  is  already  done  ,n 
the  papers  by  Hienvenu  and  Kopp  f  1 0,20,2'}. 27].  as  well  as  by  others  ' J < i 1 . 

If  the  discovery  of  a  noise  matrix  best  lilted  to  the  real  world  should 
prove  to  lie  difficult,  a  generalization  could  be  suggested  where  G  would 
be  replaced  bv  a  sum  of  weighted  noise  matrices  corresponding  to  diflereut 
st  met  ures  of  uncorrclaled  noise.  For  example 

CP  lmmI.V  CT |  G  |  ( (I  |  .  I>\  )  fTn G;(ll;,b|  1-1.7  | 

is  one  possibility  in  which  all  the  free  parameters  could  be  adjusted  to 
minimize  the  rank.  Again,  the  success  of  representing  the  noise  matrix  in 
this  way  is  reflected  by  the  degree  of  rank  reduction  in  matrix  C,.  As  we 
try  to  increase  the  number  of  noise  matrices,  a  natural  limitation  should 
appear  in  the  finding  of  insignificant  scale  factors  rr.  a  test  that  we  have 
exhausted  the  possjbjlif ies  of  noise  representation. 

Any  way  t  he  field  is  wide  open  to  improve  such  methods  or  suggest  ol  her 
ones,  and  experience  will  probably  bring  reliable  simplifications.  For  those 
who  immediately  need  any  method,  the  only  resource  is  to  carefully  study 
the  eigenvalues  of  C’P  listed  in  decreasing  order.  The  smallest  ones  may  be 
so  small  that  they  could  be  considered  negligible.  In  of  her  cases  one  can  de¬ 
tect  a  sharp  difference  between  a  subset  ol  large  eigenvalues  and  a  subset  of 
small  ones.  P  would  then  be  the  number  of  large  ones.  If  need  be,  nothing 
prevents  us  from  building  a  parameterized  model  ol  the  generalized  noise 
matrix  which  is  dependent  on  several  free  parameters,  i.e.,  rrGfu,  />.... ). 
What  has  been  suggested  requires  minimal  assumptions  because  only  the 
general  form  of  the  noise  spatial  structure  has  to  be  postulated.  Neverthe¬ 
less.  as  mentioned  in  section  d.l.  these  assumptions  have  to  lie  obtained 
from  total  ignorance  using  a  priori  knowledge. 


4.5 


The  procedure  for  determining  the  matrix.  of  the  fully  coherent  spatial 
sources  requires  a  parametric  model  to  eliminate  I  lie  second  order  noise 
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statistics.  Any  further  processing  of  matrix  C,  can  now  be  qualified  as 
"high-resolution''  with  regard  to  the  “perfect1'  sources.  Matrix  C,  is  hound 
to  the  source  vectors  S(,  by  (1.1).  The  problem  is  to  determine  whether  or 
not  w e  can  explicitly  find  the  vectors  S(>  when  C,  is  known.  See  references 
|s;_>  |  . 

Since  the  order  of  C,  is  A  and  its  rank  is  /’.  there  are  /’  nonzero  eigen¬ 
values  and  the  corresponding  eigenvectors  describe  a  /’-dimensional  sub¬ 
space  of  the  A  -dimensional  space.  VVc  will  now  show  the  wavefront  vectors, 
S,,.  also  span  I  he  same  /’-dimensional  subspace.  In  the  (  Ar  -  P )-dimensional 
complementary  subspace  the  eigenvalues  are  all  nulls  so  the  associated 
eigenvectors  are  indeterminate  which  we  shall  find  to  be  unimportant. 

The  set  of  nonzero  eigenvalues,  in  decreasing  order,  is: 

*  i  •  <  2  •  •  •  • < . . f  I’ 

and  the  corresponding  normalized  eigenvectors  (each  a  column  matrix  with 
A  components,  like  the  S;))  are: 

»ii .  do  . . .  a,,. ...  Tip 

with  the  usual  orthonormal  conditions: 


tip  ii7  -  0,  p  p  q. 

latrix  C,  may  be  written  in  canonical  form  as  justified  in  appendix 


c,  vfp  iipuj  , 


here  the  lipii,,  are  dyadics  of  the  vector  iip.  Comparing  with  (4.1)  we 


can  write 


^  SpSj  V]  <r  [uptijl  : 


i  i 


both  'ides  are  the  same  matrix.  II  we  now  look  at  the  transform  of  vector 

u  • .  we  have 


.  / 

s,.  (  S;t  ti,  )  ^  Iptip  (tijtii  )  . 


(4.10) 
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where  scalars  are  shown  between  parenthesis.  Because  of  the  above  or- 
i honormality  conditions  between  these  ii;,,  the  right  side  o(  (1.10)  is  reduced 
lo  (([pi,.  Therefore  vector  li t  appears  as  a  linear  combination  of  the  vec¬ 
tors  S,,.  As  a  consequence  li,  in  c<  asiirih/  lx  lotnjs  to  tin  suit.*  pm  ■<  .1  pun  mil 
hi i  S.,.  1  he  same  could  be  sliown  for  every  ii;).  Conversely  the  unknown 

S.,  In  Inin/  to  tin  P -ilinnn.iiomil  jtuhspacc  ihsrrihnl  lit/  the  known  r i;) .  This 
is  the  first  important  feature  of  the  vectors  S;>:  it  is  not  the  only  one.  In 
this  subspace  the  S;,  are  bound  by  several  other  relations  to  be  discussed 
next .  lncidentlv.  the  above  discussion  disproves  a  possible  misconception 
that  individual  sources  can  be  described  bv  individual  eigenvectors. 


Things  will  be  easier  if  we  accept  a  preliminary  ‘“theorem"  which  is  demon¬ 
strated  in  appendix  l'\  bet  us  consider  a  set  ol  P  vectors  V;,  describing  a 
/’-dimensional  subspace  within  an  A  dimensional  vector  space,  and  also  a 
diagonal  matrix  II  of  order  A.  built  as  follows: 


P  rows 


A'  P  rows 


ie  matrix  II  mav  now  be  written  as  a  sum  ol  dyads. 


vv  V' 
v  >'  /•“ 


(1.12) 


here  vectors  V„  have  the  following  features: 


I.  Only  the  first  P  com  p<  >  i  i  <  •  n  t  >  o|  ea<  h  \  dilhr  Irom  zero.  In  otlu'i 
words  vectors  V,,  describe  a  /’-dimensional  subspare  spanned  by  the 
first  P  axes  ol  the  lull  A -dimensional  space. 

’2.  The  v;,  will  be  <  hosen  to  satisfy  ort  honormality  comlit  ions. 
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This  being  accepted,  we  know  that  the  matrix  C,  can  be  written  as 


C,  -  UEU^ 


(4.13) 


(see  appendix  (’)•  "here  E  is  the  diagonal  matrix  of  eigenvalues,  cp,  in 
decreasing  order  with  the  last  (A’  —  P)  ones  being  null,  and  U  is  the  N- 
dimensional  unitary  matrix  whose  columns  are  the  normalized  eigenvectors 
of  C,.  set  out  in  tl  e  same  order  as  the  corresponding  eigenvalues  in  E. 
Although  the  last  (A  /’)  columns  are  indeterminate  they  are  arbitrary 

eigenvectors  and  obviously  U^U  =  I,v  a.  relation  specific  to  unitary  matri¬ 
ces.  (Iv  is  the  A -dimensional  identity  matrix.) 

Since  the  matrix  E  is  not  invertible  let  us  consider  a  diagonal  matrix  of 
real  numbers,  D  D  ,  obtained  from  E  by  replacing  the  (N  -  P)  diagonal 
zeros  of  E  with  some  arbitrary  scalar,  here  chosen  to  be  1.  It  is  clear  that 
we  can  write. 

E  -  ft  flft  (4.14) 

•  -  1 

w  1 1 ere  D  is  partly  composed  of  the  eigenvalue  square  roots.  Hence,  fol¬ 
lowing  (  1.13). 

C, -- UDMId'u^.  (4.15) 

which  may  be  solved  for  II  as  follows: 

tjtc.u  =  d-Tid'. 

D  ‘u^UEH  -  it  (4.16) 


lint  we  know  t  hat . 


C,  =  XS„Sp,  from  (4.1). 


lien  since  D.U  and  II  are  common  for  all  Sr.  we  have 
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D  'ufS 


SJUD  ' 


I  .el 


D  -  U^,. 


(4.17) 

(4. IN) 
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1  hen  we  have 


Vv„vt. 

z - ,  P  p 


(i.  in) 


Therefore  the  V,,  have  all  the  features  listed  in  items  1  and  2  at  the 
beginning  ol  this  section.  We  can  then  derive  features  of  the  S(,  by  solving 

(1.18).  Since  -  U  '1  (1.18)  yields 


S,,  :i  UI)'V(1,  (1.20) 

which  expresses  rwry  Sp  as  a  lunction  of  C,  and  of  an  A  -dimensional  vector 

V,.  gifted  with  the  features  discussed  above  in  items  1  and  2. 

Since  only  the  first  P  components  of  V;,  differ  from  zero,  the  same 

—  ,  -•  1  —  ■ .  1 
must  be  true  for  vector  Vp  D-'VP.  The  diagonal  elements  of  D-  from 

the  (P  \  I  )"'  to  the  A1’1’,  those  which  had  been  arbitrarily  chosen,  in  fact 

disappear.  Matrix  U  is  made  of  column  vectors  which  are  the  normalized 

eigenvectors  of  C3.  The  indeterminate  eigenvectors  are  neutralized  because 

when  multiplying  U  by  the  vector  Vp.  the  elements  of  U  from  the  (7’  |  I  )"’ 

to  the  Arlh  column  are  always  mult ipliecl  by  the  zero  elements  of  vector  Vp. 

Eipiafion  (  4.20)  may  be  reduced  to  its  essential  components, 


(1.21) 


vhere: 


IT  is  the  rectangular  matrix  (A'  rows  '  P  columns)  whose  columns  are 
the  first  P  normalized  eigenvectors  of  C,. 

.  .  t  i 

A  is  a  square  (  P.  P)  matrix  whose  diagonal  elements  {<■,-  ...r/,},  are  the 
square  roots  ol  the  nonzero  eigenvalues  of  arranged  in  decreasing 
older  and  corresponding  to  the  eigenvectors  of  U. 

Xp  is  a  P-component  vector  (column  matrix)  whose  components  are  the 
first  P  components  of  Vp.  The  set  of  Xp  is  gifted  with  the  same 
features  listed  in  items  1  and  2  as  the  set  of  V,,. 

Equation  (  1.21  )  is  coherent  in  the  sense  that  both  (he  left  and  right  sides 

should  be  column  matrices  with  A  rows  i.e.  dimension  (A  ■  I).  Indeed: 

Ur  is  (,V  •  PI. 

Ar  is  (/’  •  />). 

X„  is  IP-  I): 
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so  the  final  result  is  (N  x  1)  as  expected. 

There  is  a  hit  more  to  say  about  (4.21).  When  computing  the  eigen¬ 
vectors  of  a  C’SDM,  we  can  get  each  of  them  only  up  to  a  phase  factor 


exp(  iop).  If  we  choose  one  possible  determination  of  UP.  denoted  by  Up  , 

— *  p 

the  most  general  expression  of  U  will  be 


U 


r 


vP0*P 


(4.22) 


•  •  p 

re  4’  is  the  diagonal  ( P  -  P)  matrix  of  the  phase  factors.  Particularly 


pi*7  )t([>/  ff,  where  Ip  is  the  identity  matrix  in  a  P-dimensional  space. 
This  gives  (4.21  )  t he  form  S;)  -  »  L _ A  iXp 


P  p  — *  D 

u;«t>'A'xp,  or  permuting  the  two  diagonal 


mat  nres. 


Sp  -  U0P A;V'XP. 


(4.23) 


Now  the  specific  features  of  Xp  necessarily  react  on  the  vector  Zp  =  ‘p  Xp. 

1.  then  Zjzp  =  Xj  $FJ  Xp  =  1. 

t 


Since  xJX,, 


Since  Xjx,  0,  p  q.  then  zjz,  =  xj  (  ($P) T  $ 


P  v 


X,  =  0. 


lose  important  features  have  been  conveyed  from  vectors  Vp  to  vectors  Zp 
rough  vectors  Xp.  Therefore.  Zp  is  the  p^'-cohunn  of  a  ( P  x  P)  unitary 


nmtn.r  in  the  P -dimensional  subspare  of  the  first  P  eigenvectors  of  the 

rs 
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U0PAPZP 


(4.24) 


which  expresses  every  Sp  as  a  function  of  C5  and  a  partly  arbitrary  vector, 
gifted  wit  h  the  features  mentioned  above  in  items  1  and  2.  The  fundamental 
equation  (1.24)  ran  be  visualized  as  below: 
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Let  us  have: 


»z 

t  he 

»,h 

component  of 

Ha-, 

with 

1  <  n  <  /V 

~Pk 

t  he 

k'h 

component  of 

zP- 

with 

1  <  k  <  P. 

•sj»m 

t  he 

»i,h 

1 

component  of 

Sp, 

'ta 

— 

Then  (1.24)  heroines 

r 

V>  -  H  ul<rt':pk 

k-  1 


(1.2.4) 


One  ran  remark  that  the  n'1’  component  of  the  Sp  cle|)ends  on  all  the  com¬ 
ponents  of  vector  Zp.  Therefore  any  spt,  depends  on  all  arbitrary  scalars 
which  could  he  found  in  vector  Z(,  when  building  the  arbitrary  unitary 
matrix  Z.  Furthermore.  sp„  also  depends  on  all  the  eigenvalues. 


4.7 

The  first  conclusion  to  he  drawn  from  (1.24 )  is  that  an  exhaust  ive  use  of  I  lie 
('SI)M  is  not  enough  to  determine  t  he  P  source  vectors.  What  is  lacking  is 
a  (P  -  /’)  unitary  matrix  Z  whose  columns  are  the  partly  arbitrary  vectors 
Z;>  which  are  necessary  to  determine  Sp.  Using  a  priori  information  we  have 

to  choose  a  matrix  Z.  such  that  Z^Z  =  Ip. 

One  may  wonder  why  a  unitary  matrix  is  needed.  The  physics  of  this 
point  is  illustrated  by  figure  B.l  in  appendix  B.  As  mentioned  then’,  if 
♦  he  matrix  of  filters,  G.  is  given  by  the  normalized  eigenvectors  of  Cs,  the 
P  nonzero  outputs  (the  zero  ones  yield  nulls)  are  not  correlated  and  their 
spectral  densities  are  respectively  the  known  eigenvalues.  On  the  other 
hand  we  remember  that  the  P  sources  are  postulated  to  he  independent. 
That  could  lead  one  to  think  that  we  have  at  last  disentangled  the  sources 
Irom  each  other.  Unfortunately  this  is  not  trim  in  general.  A  pair  ol  uncor¬ 
related  out  puts  can  bethought  ol  as  t  wo  mixt  ures,  each  with  contributions 
Irom  all  ol  the  sources,  in  which  partial  correlations  coming  from  of  one 
source  occur  in  both  mixtures.  But  it  happens  that  these  partial  corre¬ 
lations  exactly  cancel  each  oilier,  fherelore  each  eigenvalue  represents  a 
combination  ol  all  /’  source  spectral  densities. 


t  .VA'.VO ir.v  wavefronts 
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Now  since  the  eigenvalues  are  known,  there  is  a  posibility  to  go  a  bit 
further  by  whitening  all  the  P  outputs,  thus  flattening  the  spectrum  to  a 
uniform  density.  This  is  the  purpose  of  matrix  AP  in  (4.24).  Thus,  the 
most  we  can  obtain  using  only  the  CSDM,  are  P  sped ral- normalized  mix¬ 
tures  of  the  sources.  That  explains  how  the  separation  of  sources  requires 
assumptions  which  can  be  expressed  as  a  unitary  or  “rotation”  matrix. 

One  case  where  a  solution  can  be  reached  without  additional  assump¬ 
tions  occurs  when  P  -  1,  where  obviously  the  wavefront  is  the  eigenvector 
i  i  i .  the  spectral  density  c,.  and  the  only  source  S].  The  matrix  Z  is  then 
reduced  to  its  first  element,  whose  modulus  is  equal  to  1,  and  therefore 
corresponds  to  an  arbitrary  phase  factor.  This  is  of  no  consequence  since 
a  wavefront  represents  the  relative  phases  of  its  components  and  is  not 
changed  by  a  phase  facor  common  to  all  the  elements. 

Anot her  bold  case  where  a  solution  appears  from  (4.24)  occurs  when  the 
array  is  so  large  and  super-resolving  that  all  the  source- vectors  are  virtually 
or)  hogonal.  Then  the  eigenvectors  are  the  wavefronts  and  the  corresponding 
eigenvalues  are  the  spectral  densities.  Under  these  conditions  the  matrix  Z 
is  reduced  to  a  diagonal  matrix  of  arbitrary  phase  factors,  one  per  source 
vector. 

Before  discussing  the  general  case  it  might  be  interesting  to  check  that 
I  1 .2  I )  is  com  pat  ible  with  (1.13)  when  the  useless  components  of  the  ( N  —  P  )- 
dimensional  complementary  subspace  are  removed  from  matrix  C3.  If  we 
express  S(,S(>  in  the  same  form  as  (4.24)  we  get 

s,.sj  Ui,ArZpzjA/,(U0p)t.  (4.26) 

The  matrix  Z„Z l  is  the  dyadic  of  Zp.  Now 

p  _ 

c,  l  s„s; 

;>  i 

UrAr 

I  lie  matrix  between  brackets  is  the  itlfiilili/  matrix  of  order  P.  Imleed  any 
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column  vector  v  is  transformed  as 


(!.2S) 

The  scalar  (Z^v)  is  precisely  the  "projection”  ol  v  on  the  unit  vector  ol 
the  /i'1'  axis.  So  (-1.2S)  simply  means  the  transform  of  v  is  the  vector  sum 
of  its  own  components  along  the  P  axis,  and  therefore  v  itself.  Hence 

C,  -U0r[Ap]2[U0P]f,  (I.2H) 

where  [Af  p  is  the  diagonal  matrix  of  the  eigenvalues  E  in  the  only  subspace 
where  they  are  not  nulls.  U„  is  the  unitary  matrix  U  reduced  to  (he  useful 
eigenvectors.  Now  we  note  the  form  of  (4. Id)  and  (1.29)  are  very  much  the 
same. 

4.8 

The  problem  is  now  to  correctly  assess  the  degrees  ol  freedom  required 
by  the  “arbitrary"  unitary  matrix  Z  in  order  to  determine  the  necessary 
a  priori  assumptions.  While  enough  assumptions  are  needed  to  compute 
the  Sp  using  the  full  power  of  the  array,  we  wish  to  avoid  overconstraining 
the  problem  to  preserve  the  information  contained  in  the  data  CSDM.  To 
do  this  a  unitary  matrix  ( /’  ■  /’)  is  built  of  P2  complex  scalars  which  are 
equivalent  to  'IP2  real  scalars.  The  constraints  between  these  scalars  are: 

1.  Every  column  must  be  a  normalized  vector  giving  one  real  scalar 
relation  per  Sp  for  a  total  of  P  relations. 

2.  Every  column  must  bo  ort  hogonal  to  all  the  ot  hers,  which  corresponds 
to  /’(  P  1  ) / 2  relations  between  complex  numbers  and  therefore  a 
total  of  /’(  P  1  )  relations  between  real  scalars. 

So  the  2P:  real  scalars  ol  Z  are  bound  by  P!  relations,  and  the  arbitrary 
freedom  within  matrix  Z  would  seem  to  be  dependent  on  /’”  rial  scalars. 
However  we  have  seen  that  the  phase-fac  tor  <il  each  column  Zp  of  Z  plays  no 
role  in  our  analvsis  thus  allowing  <>in  real  scalar  to  remain  undetermined. 
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As  a  result  the  number  of  real  scalars  needed  to  get  all  the  Sp  is  P2  —  P  = 
/’(/’  1). 

If  we  look  at  the  set  of  P  vectors  given  by  (4.24),  each  corresponds  to  N 
complex  scalar  equations,  or  equivalently  to  2 N  real  scalar  equations.  For 
the  whole  set  we  have,  2N P  real  scalar  equations.  What  is  the  number  of 
the  unknown  parameters?  We  already  know  there  are  P(  P  1)  real  scalars 
introduced  by  the  vectors  Zp,  gathered  in  matrix  Z.  On  the  left  side,  all 
2.V/’  scalars  forming  the  components  of  the  Sp  are  unknown. 

We  then  have  more  unknown  parameters  than  equations.  That  is  why 
Ini  ill  ignorance  about  the  Sp  prevents  any  solution  to  emerge,  even  with 
the  most  extensive  use  of  the  C-SDM.  Something  has  to  be  assumed  about 
the  SP:  but  any  assumption  should  provide  precisely  P{  P  —  1)  real  scalar 
iijiialiims  between  the  components  of  the  Sp.  More  assumptions  would  lead 
us  to  impose  specific  features  (relations  between  the  elements)  on  both  the 
(  SI).M  and  the  related  matrix  U0  A  ,  therefore  causing  conflicts  between 
the  experimental  results  and  the  assumptions. 

before  trying  a  general  approach  using  the  minimum  necessary  assump¬ 
tions.  it  is  interesting  to  illustrate  the  previous  development  with  the  case 
I'  2.  A  unitary  matrix  of  order  2  is  made  originally  of  8  real  scalars.  But 

the  constraints  Z^Z  I?  lead  to  the  classical  writing, 

2  (  cost,!.')  si ii ( V )exp(  id)  \  f  exp (uf>x)  0 

\  si n ( v  )exp(  id)  cos(r)  /  \  0  exp(i(/>2) 

where  only  I  real  scalars  remain.  The  diagonal  matrix  on  the  right  sirle 
corresponds  to  t  ho  phase  factors  of  the  lirst  and  second  columns.  They  can 
remain  arbitrary  so  that  only  2  real  scalars  w  and  ft  have  to  be  introduced 
m  t  he  vector  equations. 

It  must  be  observed  that  both  scalars  are  present,  in  every  element  of 
the  matrix  and  particularly  in  all  columns.  The  generalization  of  this  for 
larger  values  of  P  is  a  harbinger  of  computing  difficulties  which  arise  when 
detei  mining  the  source  vectors  as  functions  of  the  arbitrary  scalars.  In  the 
general  case  we  are  confronted  with  finding  the  number  ol  scalar  relations 
between  the  components  ol  the  Sp  which  are  necessary  to  complete  the  set 
of  equations  (  1.24)  with  exactly  as  many  unknown  parameters  as  there  are 
available  equations.  Not  violating  this  condition  with  reliable,  compatible 
assumptions  is  hopefully  a  way  to  solve  the  problem  ol  finding  source  veef ors 
as/i’i/  tht  full  anai/  pain  r. 
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For  the  case  P  —  2,  (4.24)  becomes: 


cos(v’)  \ 

si u ( (/•  )exp(  id)  ) 

-sin(  i')exp(  i9) 
cos(  </- ) 


(4. 


These  two  vector  equations  lie!  ween  the  components  of  S[  and  S2.  yield  2.Y 
complex  scalar  equations  or  4/Y  real  scalar  equations.  On  the  left  side  of 
(4.30)  the  unknown  parameters  to  he  computed  are  precisely  the  4A  real 
scalars  of  the  components.  On  the  right  side,  the  unknowns  (unavoidable 
whether  interesting  or  not)  are  ir  and  9.  There  are  two  more  unknowns 
than  relations  so  we  must  either: 


1.  Arbitrarily  forget  about  p  and  9  by  deciding  <r  0  and  assume  the 
two  source  vectors  ate  orthogonal  (very  large  arrays),  a  rather  crude 
procedure. 

2.  “Invent"  2  real  scalar  relations  between  the  2 A  real  scalars  of  the 
source  vectors.  This  is  not  very  stringent  and  is  exactly  what  is 
required  to  reach  a  solution.  More  assumptions  would  lie  too  much 
and  would  require  compatibility  with  the  experimental  ('SI)M. 


4.9 


Let  us  first  remember  that  we  are  considering  arrays  of  some  (it)  sensors  or 
more,  and  P  up  to  20  or  30.  Clearly,  all  the  assumptions  made  in  previous 
sections  about  the  source  vectors  are  of  the  “overinformation  ”  type,  as 
compared  to  the  minimum  sufficient  assumptions.  For  example  t  he  “plane- 
wave  assumption”  of  chapter  2  corresponds  to  a  source  vector  structure 
with  a  complex  amplitude  and  two  angles  for  the  direction.  A  total  of  I 
real  parameters  are  needed  for  this  vector  description.  It  actually  implies 
2  Y  4  implicit  relations  between  t  be  2.Y  real  scalars  of  the  components 
giving  a  total  of  P(2N  I )  relations  lor  all  the  source  vectors  which  exceeds 
the  /’(  P  1)  strictly  necessary  relations  compatible  with  t  lie  experiment  al 
('Sl)M  C,.  Sometimes  the  minimum  relations  are  praclicallv  anticipated. 
As  an  example,  for  a  linear  arrav  <*l  equispaced  identical  sensors,  the  source 
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vector  ('SI)M  is  <>1  the  Toeplitz  type.1  The  Toeplitz  criterion  helps  separate 
i  lie  source  CSDM  C,  from  the  data  CSDM  Cr,  which  contains  unwanted 
noise,  thus  giving  an  estimate  of  P  (the  problem  is  to  lind  the  largest 
possible  Toeplitz  matrix  contained  in  C,.  ). 

Vezzosi  [25|  has  given  a  lull  solution  under  less  stringent  assumptions. 
For  a  source  vector,  only  the  amplitude  of  the  components  is  supposed  to 
be  constant.  T  his  common-value  amplitude  and  the  phases  must  then  be 
determined.  The  assumption  of  a  constant  amplitude  for  N  components 
corresponds  to  A  I  real  scalar  relations  for  each  source  vector  and  P(N -l ) 
relations  for  all  sources.  Although  this  is  a  big  step  forward,  it  is  still 
overconst  rained  since  if  exceeds  the  minimum  P(  P  —  1)  relations.  Again, 
the  source  CSDM  is  somewhat  specific  thus  giving  a.  clue  for  the  estimation 
of  P. 

Other  examples  of  this  kind  could  be  proposed  more  or  less  based  on 
intuition  about  specific  situations.  True,  it  is  painful  to  refrain  from  sim¬ 
plifying  an  enormously  diflicult  problem  by  adjusting  the  CSDM  to  a  given 
shape,  and  then  further  refrain  from  blaming  any  discrepancies  with  the 
experimental  CSDM  upon  a  lack  of  spectral  and  cross-spectral  density  es¬ 
timation  accuracy.  CSDM  shape  adjustment  becomes  even  more  tempting 
because  it  yields  an  estimate  of  P,  and  we  must  admit  that  in  section 
l.l  we  could  only  derive  such  estimates  from  conjectures  about  the  set  of 
eigenvalues.  But  this  is  only  transferring  conjectures  from  one  place  to 
another.  Clearly  compromises  may  be  successful  in  some  situations.  Still, 
adjusting  a  problem  to  its  solutions  leaves  us  uneasy.  Anyway,  for  better 
or  worse  we  have  chosen  to  sail,  so  we  shall  now  try  to  lind  a  minimum  set 
ol  assumptions  which  balance  the  equations  with  the  unknowns. 


t:' 


4.10 

( Ini  goal  is  now  to  suggest  something  about  the  source  vectors  or  the  wave- 
fronts  themselves  (and  not  about  a  mode)  of  the  medium  since  we  do  not 
pietend  to  localize)-  It  should  be  a  “minimum’  assumption  which  is  both 
complete  and  reliable  thus  allowing  us  to  use  the  full  power  ol  the  array  as 

1  A  o.H-lriint  impos.<|  on  tin-  C.  matrix  lw  I  lie  "plane- wave  assumption"  made  for  any 
V  senior  arrav.  lias  some  internal  order  snrli  as  tin  loeplitz  matrix,  nnlv  more  < I i (lie r 1 1 1 
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expressed  in  ihe  (’Sl)M  through  the  set  of  vector  equations  (4.21). 

l'lie  assumption  proposed  here  is:  a  source  vector,  which  is  not  nec¬ 
essarily  n  plant  trace,  can  he  dt  scribed  as  a  coherent  sunt  of  plane  warts 
across  tht  a  pert  art  of  the  array.  Such  a  superposition  should  approxi¬ 
mate  the  source  vector  more  closely  when  the  number  of  plane  waves  Is 
large.  hut  this  number  is  limited  by  the  number  ol  sensors  in  the  array.  A  . 
This  assumption  may  he  disputed,  but  it  seems  to  deserve  the  qualification 
"minimum'  '  be  cause  as  we  s  hall  see  later: 

1.  The  number  of  plane  waves  is  act  idly  large  in  most  of  the  situations 
we  can  consider. 

2.  The  geometrical  significance  of  a  plane  wave  is  not  here  an  absolute 
preliminary  requirement  but  is  only  to  be  a  posteriori  controlled. 

The  last  remark  gives  the  wavefront  model  large  flexibility.  If  an  ex¬ 
pected  "plane"  wave  is  found  to  he  non-plane  (we  shall  see  how),  it  still 
contributes  to  a  parametric  model  of  the  wavefront.  If  the  plane-wave  fea¬ 
ture  is  confirmed,  it  gives  a  clue  for  localizat  ion.  That  is  why  we  shall  name 
the  model  a  coherent  sum  of  “conditional  plane  waves."  CPVY. 

flunking  of  large  fixed  arrays,  let  us  now  restrict  ourselves  to  an  array 
with  .V  point-like  sensors  (as  compared  to  t he  smallest  wavelength  involved) 
which  have  a  flat  frequency  response  over  the  useful  bandwidth.  Since  the 
array  geometry  is  completely  known  we  can  define  a  reference  system  in  the 
real  three-dimensional  space  with  a  point  origin  and  three  orthogonal  axes. 

I  he  Cartesian  position  of  the  ;?11'  sensor  is  then  given  by  a  known  vector 
which  specifies  the  length  and  direction  from  the  origin  to  the  sensor.  f„. 

On  the  other  hand  an  incoming  plane  wave  is  completely  described  by 
a  normalized  3-dimensional  vector  n,  |o’|  -  1,  pointed  from  the  origin 
toicards  the  plane  wave  (opposite  to  its  propagation  direction)  representing 
two  angles  (bearing  and  tilt).  This  vector  tlocs  not  depend  on  frequency 
and  is  unknown  in  our  problem.  A  “complex  amplitude’’  .4  representing 
the  amplitude  and  phase  at  t he  origin  point  is  also  unknown  and  frequency 
dependent  . 

\  cry  classically  the  phase  factor  corresponding  to  the  travel  time  delay 
from  the  origin  to  the  sensor,  is  exp  2rrda.\,  •  o  )j  where  the  dot  stands 
lor  a  scalar  inner- product  and  where  the  vector  is  the  vector  r,,  men 
s u red  m  terms  o|  wavelengths.  I  lie  "length  >4  vector  v,,  is  proportional 
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i<>  frequency  but  remains  completely  known  in  every  frequency  bin.  As  we 
have  ilone  throughout  this  paper,  we  let  the  frequency  dependence  be  im¬ 
plicit.  Therefore  is  a  known  vector  in  the  d-dimensional  space,  while  for 
a  "true  plane  wave,  o  is  unknown,  frequency  independent,  and  represents 
two  real  scalars,  For  a  CPW.  vector  d  will  be  computed  in  every  frequency 
bin.  Only  a  posteriori  will  the  CPW  be  verified  as  truly  being  a  frequency 
independent  plane  wave. 

As  for  the  complex  amplitude  .1,  obviously  |.1|  is  the  same  on  all  sensors 
and  at  the  origin.  If  we  were  to  represent  only  one  plane  wave,  .1  could  be 
a  real  positive'  scalar,  and  the  above  phase  factor  wotdd  express  the  phases, 
relative  to  that  of  the  origin.  But  since  we  shall  deal  with  a  sum  of  coherent 
plane  wave's,  we  have  te>  describe  the  relative  phases  of  several  plane  waves 
at  the  origin,  anel  .1  has  to  be  a  complex.  frequency  dependent,  scalar.  Se> 
then  1  ile’pe'tiels  on  2  real  scalars,  amplitude  and  phase,  both  unknown  in 
e •  v •  •  r \  frequency  bin.  I'he  expre'ssieui  for  a  plane  wave  on  the  n1'1  sensor  is 
t  lie'll 

. lexp  [2 ni  ( u.'„  •  d  )) 

and  depends  on  I  real  scalars.  While  two  ol  them  .should  not  be  frequency 
•  h'peiuh'iit .  this  may  not  always  occur. 
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Auv  'onree  vector  S  described  as  a  sum  ol  A,  coherent  plane  waves,  yields 
on  the  o'1'  sensor  a  component. 

h . 

^  Ak exp  [2 ni  (u3n  •  )]  . 

h  -  1 

For  some  particular  source  vector  S;,  with  a  specific  number  of  CPW,  Kp, 
tin-  n'!'  component  is  given  by 

K, 

S,„,  ^  .C,,exp  '2t;  ( u\,  ■  or(,  )i  :  .  (1251) 

!■  i  1 

l\ ..  ma\  not  be 1  tin'  same  b>r  all  Sp  In-cause  <-a<  h  source  may  have-  a  diflerent 
in  nil  i  pa  1 1 1  a  rri  va  1  si 1  int  ii  re.  1  Ins  album  u^  flexibility  in  the  way  we  express 
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anv  loose  a  priori  idea  we  might  have  about  the  relative  tomplexity  o|  the 


source  vectors. 


component  .s  ,  depends  on  \  l\ ,,  r 


eat  scalars. 


we  follow  (1.21  )  from  one  sensor  to  another  only  the  geometric  vector  u.'\, 
varies,  hut  it  does  so  in  a  perfectly  known  way.  1  herelore  for  each  ol  1 1 1« •  S., 
components,  we  have  a  set  of  \l\p  unknown  real  scalars.  The  phase  of  .  h( 
for  a  given  k .  is  the  same  on  all  sensors  and  expresses  the  relative  phase  of 
the  k,u  (T\Y  huilding  S,„ 

Since  an  S,,  cannot  he  determined  better  than  within  a  phase  factor,  we 
could,  in  principle  take  an  arbitrary  phase  for  one  CI’W  .  for  example  I, 
This  would  reduce  by  one  the  number  ol  unknown  real  scalais.  However  we 
shall  not  do  so.  because  we  are  going  to  compare  t  his  paramet  ric  expression 
of  S;>  with  that  given  by  (  1.2-1)  or  (1.2').  These  two  expressions  must  only 
differ  bv  an  unknown  phase  factor  since  they  each  represent  the  same  S., 
except  for  the  phase.  Thus  we  must  maintain  a  free  phase  factor  lor  l|(,  to 
represent  anv  possible  phase  difference  between  the  two  expressions  for  the 
same  S,,.  As  a  result  the  number  of  free  real  scalars  in  the  CI'W  model  is 
confirmed  at  M\r.  The  total  of  CPW  for  all  the  S,,  is  l\  with 


l'  /V 


(  1.221 


and  tli<‘  problem  is  to  assess  the  value  of  /\  by  coni  parison  wit  li  (  1 .2  I  I . 
latter,  as  w<'  remember,  alreadv  involves  /’ i /’  I  )  free  real  scalars. 


4.12 

At  the  end  of  section  l.<>.  we  have  already  noticed  that  the  n’1  compoiu'iit 
of  Sp  from  (1.21 )  in v< lives  all  the  component s  of  Z(1  and  t here lore.  the  same 
/’(  /’  I)  free  real  scalars  of  matrix  Z  (obvious  in  the  /’  2  example  of 
section  !.*).  Now  the  comparison,  for  the  same  S;,  ol  the  ('SI)M  derived 
expression  and  the  ('I’W  model  yields: 

I  vert  or  eipia  t  ii  >n 
A  i  omplex  scalar  eipiat  ions 
2 A  real  scalar  eipiations. 

for  the  whole  set  of  S;,  we  have  2 A  /’  real  scalar  eipiations.  I  lie  number  ol 
free  real  scalars  to  be  derived  front  this  sei  o|  eipiations  is  |  /\  .  /' |  /’  I  i 
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nml  lln'  value  of  A  is  given  by  I A  I-  P{  P  1  )  2  AT,  or  solving  lor  A  , 


XP  P(P  I) 


(4.33) 


I  lierelore  lor  a  given  /’.  1  lie  number  of  available  C’PW  is  obviously  growing 
uilh  A  (  (he  array  power)  at  the  rate  of  P  waves  for  two  more  sensors, 
hi  me,  one  more  wave  per  wavefront. 

Perhaps  less  obvious.  A  is  also  ijrounmj  with  ihe  number  of  sources  /’, 
within  the  limits  of  its  possible  variation  for  a  any  given  A  .  Indeed  the  P 
derivative  of  the  right  side  of  (1.23)  is  (A’  -  P)i'l  +  I /l,  a  positive  definite 
number.  It  is  an  important  aspect  of  this  assessment  that  the  whole  set  of 
Miiitee  vectors  is  described  better  when  there  are  more  sources.  Later,  we 
shall  come  back  to  this  remark.  Nevertheless  the  mean  value  of  OPW  per 
.oM/ni  decreases  as  P  grows.  Though  we  have  mentioned  that  the  CPW 
mas-  not  be  land  cannot  always  be)  equally  distributed,  a  rough  estimate 
* >1  the  power  of  the  array  per  source  vector  is  given  by  the  integer  nearest 


(•1.34) 


(  «  » 

A'  A’  P  -  1 

P  2  "  T 

Tim  number  decreases  with  P.  and  the  P  derivative  of  the  right  side  is 
I  I  wlibh  means  that  for  1  move  sourer  s  t  lie  description  of  each  individual 
wavefront  is  reduced  hi /  only  one  CPW.  an  exremely  slow  rate. 

[fefore  giving  examples,  we  shall  first  make  aim!  her  small  assumption; 
A  r-  ossunntl  to  in  mi  inn  number.  Considering  arrays  ol  GO  sensors  or 
more  this  is  anything  but  stringent.  Equation  (1.3.3)  reads 


A 


(4.3  b) 


with  A  2  an  integer,  for  A  to  also  be  an  integer,  the  even  number  P(P  -  1  ) 
must  be  a  multiple  of  1. 

I  et  us  see  a  few  examples  where  this  condition  is  fulfilled  and  afterwards 
hat  happens  when  it  is  not : 


I  A  lit).  /’  20  hence  A  b()b 

I  he  whole  set  of  wavefront  s  is  represent  ed  by  btlb  ( '  P\\  .  I  here  ca n  n<  >1 
be  the  same  number  of  wavefronts  per  source,  but  the  mean  power  ol 
t  In-  a  n  a  v  <an  be  represent  ed  b\-  the  integer  n  ea  rest  t  o  btlb  20.  namely 


t! 
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2.  A  -■  t>2.  P  —  20  hence  A  —  5 20 

Adding  2  more  sensors  (A  an  even  nnmher),  the  whole  set  of  wave- 
fronts  are  now  described  l>y  20  more  CPW.  one  per  source.  Indeed 
t lie  integer  nearest  to  525  20  is  20. 

3.  A'  -  00,  P  21  hence  A  -  532 

Coming  back  to  00  sensors  but  with  4  more  sources  as  compared 
with  case  1,  the  global  situation  is  described  by  77  more  C|’\V.  More 
sources  give  more  information  as  could  be  expected.  Nevertheless, 
per  source,  the  integer  nearest  to  532  21  is  now  21.  (I  is  remarkable 
that  4  more  sources  have  stolen  only  one  CPW.  as  compared  to  case 
1.  for  the  descripton  ol  a  particular  source  vector. 

4.  A  -  00.  P  22  hence  A  -  545  1  2 

Let  us  take  now  an  intermediate  example  where  the  even  number 
/’(  P  1)  is  not  a  multiple  of  1.  For  this  case  the  integer  nearest  to 
545  22  is  either  24  or  25  and  the  (lower  of  an  individual  wavefront 
description  is  about  the  same  as  in  cases  1  and  3. 

What  is  the  meaning  of  a  "half-CPW"?  It  means  that  515  CPW  is  a  bit 
too  much  and  actually  brings  into  the  system  of  real  scalar  equations  too 
many  free  real  scalars  for  the  number  of  equations.  One  of  the  5  15  CPW 
has  become  "half-known. "  depending  on  only  2  free  real  scalars  instead  of 
4.  li  r  then  nrnl  some  additional  a  priori  informal  ion .  One  can  build  om 
of  the  545  CPW  with  a  known  direction  vector,  o',  which  is  not  frequency 
dependent,  a  tru<  plane  wave.  The  complex  amplitude  remains  undefined, 
lo  choose  the  direction  we  can  rely  on  any  loose  a  priori  idea  we  might  have 
on  the  significant  presence  ol  a  particular  piane  wave  in  the  description  of 
one  of  the  wavefronts.  It  is  not  very  critical.  An  error  would  be  revealed  by 
a  small  corresponding  amplitude  in  the  final  computed  solution.  One  could 
say  that  the  stringent  assumption  of  chapter  2  has  now  been  submerged  as 
a  very  small  part  within  a  flexible  solution. 
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The  number  of  free  real  scalars  available  to  describe  one  wavefront  is 

I  A 

2 A  P  ■  \.  i  1.3b) 
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It  does  not  seem  preposterous  to  say  this  number  is  “large.”  In  the  previous 
examnles  it  is  about  100,  more  precisely  it  is: 

101  in  case  I . 

I  O')  in  .  ase  2. 

07  in  case  0, 

00  in  case  1. 

This  sort  of  assessment  does  not  help  us  at  all  in  determining  the  clustering 
o|’  tlmse  scalars  for  a  sequence  ol  coherent  CI’VV.  But  it  does  help  us  obtain 
a  rough  evaluation  of  the  parametric  processing  power.  For  an  efficient  GO 
sensor  array,  it  seems  sensible  to  associate  a  number  of  “perfect”  sources 
between  10  and  10.  l  lie  comparable  figures  are  then: 

1  I  I  free  scalars  for  It)  sources, 
s  I  free  scalars  for  If)  sources. 

which  exhibit  a  ratio  of  only  0.73  scalars  for  I  times  more  sources.  Is  it. 
going  too  far  to  say  there  is  apparently  never  a  shortage  of  free  scalars  and 
their  number  varies  only  slowly  with  the  number  of  sources?  Both  features 
are  those  of  a  good  parametric  representat ion. 

.Nevertheless  it  must  be  emphasized  that  in  order  to  process  a  set  of 
|l!  ounces  according  to  the  procedure  in  section  1.1.  we  must  first  find 
a  parametric  noise  model  (no  wavefront)  able  to  “eat"  5(1  eigenvalues  out 
of  the  fit)  x  GO  data  matrix  C,.  and  then  get  the  source  matrix  C,.  Of 
course  a  large  array  might  be  kind  enough  to  make  the  subspace  division 
nearly  obvious  depending  on  the  eigenvalue  spread.  Still  the  determination 
ol  /’is  the  weak  point  of  the  method,  even  though  we  might  have  some 
due,  I  ( lou hi  fill).  All  of  these  difficulties  arise  from  trying  to  avoid  faith  in 
stringent  wavefront  assumptions.  But  the  assumptions  we  have  made  are 
n  la  I  in  It/  wtak.  Instead  of  direct  statements  about  the  Sp  themselves,  we 
make  indirect  statements  about  eigenvalues,  which  are  each  only  a.  small 
part  of  some  S;,  component  [see  equation  ( 4 . 2 -r> ) ] .  Also  we  have  not  vet 
used  either  the  frequency  dimension  or  any  experimental  data  about  the 
spreading  of  the  eigenvalues  obtained  under  controlled  conditions. 

The  analysis  of  the  data  matrix  C,.  must  be  performed  in  hundreds  of 
It  e<  p  i  <  ■  1 1  c  v  bins,  (’ertainly  in  an  operational  situation  /’  is  not  expected 
to  be  the  same  everywhere  within  t  he  whole  bandwidth.  Ifitlerent  sources 
ma’.  Iiausmit  component*-  itt  various  dith-rent  sub  bands.  Nevertheless, 
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when  following  the  same  eigenvalue  across  a  lew  neighboring  liins.  some 
stability  or  smooth  variation  shouhl  he  observed.  especiallv  lor  the  large 
bandwidth  sources  considered  hen*,  hast,  erratic  lluct  nations  should  dis 
qualify  an  eigenvalue  as  being  represental  ive  of  the  set  of  S  ,. 

On  the  other  hand,  experiments  with  P  controlled  sources  are  needed 
land  are  already  being  done)  to  gel  some  data  about  the  actual  spreading 
of  the  eigenvalues  and  the  reliability  ol  various  methods  of  estimatin';  the 
value  ol  /’  (as  il  it  were  unknown).  As  small  as  we  have  left  the  part  of 
estimating  /’.  it  must  be  backed  by  experimental  data  even  though  there  is 
no  universal  and  definitive  criterion  to  be  expected.  In  short,  the  behavior 
o|  an  eigenvalue  across  a  low  adjacent  frequency  bins  should  contribute  to 
its  accept  a  nee  as  I  >eing  represent  at  i  ve  ol  a  set  of  s<  uirce  veet  ors  and  I  lien ‘lore 
indicative  ol  the  number  of  sources.  This  is  apparently  the  corner  stone  of 
t  he  eigenvalue  approach. 

I  he  a  posteriori  verification  of  a  ('P\\  as  a  true  plane  wave  may  be 
tested  by  the  stabilit  y  ol  its  direct  ion  vector  o  from  one  Irequencv  bin  to 
another.  Probably  only  a  lew  veet  ors  will  qualify  lor  t  rue  plane  -  wave  st  at  us. 
Hut  none  of  this  invalidates  the  validity  ol  the  parametric  representation. 
In  lad.  true  plane  waves  are  possibly  <>l  interest  only  lor  a  fast,  obvious 
solution  to  the  localization  problem.  Again.  across  the  bandwidth,  we 
should  observe  no  erratic  variations  of  the  a,,;  and  bound  to  source 

vector  S.,.  which  are  liable  to  call  into  question  either  the  Sr  or  the  number 
o|  source  vectors.  Indeed.  <  r<  ri/llinx/  r<li<s  on  fix  (.Pinxilxm  of  P. 

A  huge  computational  ellort.  probably  not  now  available  in  real  time, 
is  the  <  ost  to  be  paid  lor  the  implementation  o|  this  process.  I  he  basic 
mat  hemal ical  problem  to  be  solved  is  governed  by  (  l.'JI)  or  l  l.'Jo).  and 
I  l.dl  ).  expressing  the  source  vector  S;,  in  two  dillerent  ways.  \\  hen  solved 
il  yields  the  source  vector  S,.  along  with  the  apparent  spectral  density. 

•/,,  S  .S,,.  and  the  wavelront  F(,  if,,  Sp.  So  in  principle  the  problem 
ol  source  identification  with  unknown  wavefronts  is  solved  within  the  sense 
given  to  the  term  "identification  at  the  beginning  of  ihts  papei.  |.H  ns 
now  i  ••member  that  since  the  process  is  adaptive,  the  whole  solnt  ion  o  valid 
lor  a  Mine  duration  /  and  liable  to  be  refreshed  every  /  seconds  or  sooner. 
As  in  any  rich  parametric  representation  with  such  enormous  ||e\ibihlv.  the 
geometric  and  physical  signilicance  o|  so  mam  tree  paianaler-.  i-  not  v  rv 
iiii  port  a  nt . 
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4.14 

We  could  stop  this  presentation  here  since  we  are  not  primarily  interested 
in  localization.  Nevertheless  a  few  remarks  should  he  made  concerning  the 
relationship  between  what  we  have  done  above  and  the  process  of  local¬ 
ization.  First,  whatever  localization  method  is  to  he  chosen,  it  is  good  to 
have  a  line  description  of  the  source  vectors.  Second,  localization  is  not 
separable  from  some  description  of  the  propagation  medium.  In  fact  the 
source  vectors  themselves  on  a  description  of  the  medium  since  each  vector 
,  I,  ju  nils  on  Ih <  saint  mt  ilium  model  ns  well  as  the  coordinates  and  spectral 
densitv  of  the  source.  We  now  have  the  source  vectors  which  we  have  been 
li  ving  to  describe  while  using  a  minimum  of  assumptions  consistent  with 

I  lie  array  size. 

Our  present  goal  is  to  use  the  same  guide-lines  lor  localization.  We  seek 
free  parameters  which  describe  the  source  vectors  taking  into  account  the 
fad  that  the  source  vectors  are  known  functions  of  frequency.  We  shall 
approach  this  problem  only  in  a  fairly  superficial  manner.  It  must  be  made 
dear  that  what  follows  is  neither  a  method  nor  even  a  theory.  Instead  it 
is  a  conceptual  assessment  of  ( he  free  parameters  liable  to  be  involved  in  a 
parametric  description  of  the  medium. 

We  have  already  assumed  t  hat  all  sensors  exhibit  the  same  flat  frequency 
response  and  are  approximately  point-like  with  a  known  postion  given  by 
a  . !  - 1 )  vector,  r’„  representing  the  three  coordinates.  We  now  additionally 
assume  then'  is  no  calibration  problem  and  therefore,  the  symbol  S,,  which 
has  I  urn  used  to  represent  a  source  vector  in  terms  of  sensor  output  elecric 
vi  >1 1  use,  can  also  be  used  to  represent  the  sound  pressure. 

To  include  the  principle  of  reciprocity,  the  sound  field  of  a  normalized 
amplitude  source,  is  given  by  a  Green's  function 

G»  ;m  Gr  (  •!  p ,  Up .  ,  1  „  ,  Cl  i  ,  .  .  •  Cl  ^  ,  .  .  .  (I  Q  )  , 

which  is  symetric  with  respect  to  the  source  coordinates  ,rp,  f/p,  zp  and  the 

coordinates  f„  of  the  n'  '  sensor:  the  parameters  ii| - n,, _ ny  represent 

a  model  of  the  medium.  Building  the  Green's  function  is  generally  easy  lor 
vi •  i  v  simple  conditions  such  as  a  perfectly  reflecting  surface  or  a  bottom 

I I  •  1 1  ec  t  ion  with  a  known  impedance.  1  he  parameters  n,(  are  then  very  few 
and  known  a  priori.  But  the  Green's  function  rapidly  becomes  very  complex 
for  i  he  more  realistic  fields  associated  with  specific  propagation  conditions. 


(iii  i  \k\ow.x  \\'Avr:i- homs 

e\g.  <liHrac1  ion  phenemn’na  etc.  So  usually  the  parameters  and  the  hum  ol 
the  Cireeu's  function  are  obtained  Irom  an  accepted  medium  model.  Most 
of  the  t  ime  such  models  are  over  simplified  represent  at  ions  of  the  real  world 
which  only  describe  the  basic  outline  of  an  actual  acoustic  field . 

The  point  of  view  here  is  a  bit  diflerenl.  We  conjecture  that  it  might  b” 
possible  to  build  a  general  Cireeu's  function  which  involves  unknown  free 
parameters.  We  then  try  to  determine  how  many  parameters  are  needed 
for  compatibility  with  the  array  size  and  the  recently  acquired  knowledge 
of  the  source  vectors.  1  his  is  important  because  this  number  ol  available 
para  met  ers,  Q.  determines  t  he  complexity  of  t  he  Green's  fund  ion  which  can 
be  built,  l’he  challenge  facing  model  makers  and  fast  real  time  computer 
specialists  is  the  actual  building  of  such  flexible  (freon's  functions. 

Parameters  e,(  may  be  real  or  complex,  frequency  dependent  or  not .  But 
w  ha  fever  they  might  be.  t  hey  are  the  same  for  all  sources  and  sensors.  Only 
the  source'  coordinates  and  the  sensor  coordinates  vary  with  source  and, or 
sensor  posit  ion  with  t  he  la  t  ter  presumably  known.  I  hen  the  />* 1 1  source’  with 
;m  unknown  eeunplex  freupiency  elepeneh'iit  amplitude  A,,  yie-hls  on  the>  n"1 
se’ii se >r .  an  acoustic  fie-lel  A(,G,„,  which  is  precisely  the'  n11'  component  of 
ve’e  tor  S(,  now  known  from  (  f.2d).  \\’e>  must  remember  that  S(,  can  be 

known  only  up  to  an  arbitrary  phase’  laden-  which  is  that  <>l  A,,.  A I  se  > 
re’membe*r.  the-  Irnt  spectral  density  of  the’  source’  mentioned  in  chapteu- 

is  <lt  ‘A  ’’  while’  <1  S^S  was  enily  the-  apparent  spe-dral  ele’itsily  as  se ’e-n 
through  the-  array,  for  a  given  S,.  we-  have-  ,V  complex  sealar  r<’lations  of 
I  In’  type 

A  i.Tl  n  II  S 1 1  n  .  (boll 


this  redatioti  mvolvi’s  I  parameters  lienimi  to  the’  source:  1  p, t/j,.  :(1. 

which  happen  to  be’  real  scalars,  and  Q  parameters  hound  te>  the-  medium. 
(  onsieleriug  the-  whole’  set  of  /’  seturce  veefens.  we-  have-: 


rameters  beumel  to  the’  sentree 


A  I’  se  alar  relat iems. 

I/’  paramete-rs  lor  1 1 i <■  senirevs. 
<J  paramete-rs  her  the-  meelium. 

h<’  systeun  is  balanee-el  whe  n: 

I/’  •  <> 
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So.  t  lie  number  of  parameters  required  by  the  general  Green’s  function 
of  the  medium  grows  with  A  and  P  as  expected.  Nevertheless  the  full 
proportionality  to  P  in  (1.38)  is  somewhat  misleading.  It  comes  from  the 
fact  we  are  supposed  to  have  first  solved  for  the  source  vectors,  but  these 
vectors  are  already  an  implicit  description  of  the  medium.  Actually,  the 
number  Q  grows  with  P  at  about  the  same  rate  as  the  number  of  C'PW 
in  section  1.12.  So  we  can  say  that  more  sources  shed  more  "light"  on  the 
nu  ilium  description  as  could  be  reasonably  expected.  All  of  this  suggests 
tin  list  of  sources  for  stodging  the  medium. 

from  (4.38),  the  case  A’  -  60,  P  —  20  yields  Q  —  1120  which  may 
seem  an  enormous  descriptive  capacity.  But  we  must  remember  that  we 
deal  with  seniors  when  describing  complex  propagation  phenomena.  For  the 
localization  problem  we  may  only  be  interested  in  the  source  parameters.  So 
we  can  "eliminate"  the  medium  description  from  (4.37)  while  solving  for  the 
source  coordinates  and  spectral  density.  There  is  no  formal  guarantee  the 
coordinates  will  l»e  perfectly  constant  from  one  frequency  bin  to  another. 
Only  the  ultimate  fitness  of  the  model  to  the  real  situation  might  make 
it  so  as  a  first  approximation.  But  this  is  also  the  case  for  any  other 
passive  localization  process.  The  question  of  whether  or  not  this  process  is 
accurate  to  I  he  limits  of  its  high  degree  of  complexity,  is  a  mat  ter  of  further 
experiments. 

Opposite  to  the  localization  problem,  we  have  the  experimental  study 
of  tin  nxdium  with  sources  controlled  in  both  position  and  spectral  den¬ 
sity.  Individualization  of  source  vectors  may  result  from  transmissions  in 
nonoverlapping  sub-bands  or  from  the  disentangling  procedure  suggested 
in  sections  1.1  I  1.13.  Once  the  source  vectors  are  obtained,  NP  free  pa¬ 
rameters  are  available  to  describe  the  model.  Reiterating  a  remark  already 
made  about  the  continuous  plane  waves  and  large  parametric  descriptions, 
the  physical  meaning  of  a  particular  parameter  may  he  less  critial  when 
there  are  more  of  them.  Of  course,  much  further  study  will  be  necessary  to 
confirm  I  his  possibility. 


4.15 

As  a  conclusion  lor  chapter  1.  let  us  say  we  have  at  last  discarded  any  a  pri¬ 
ori  knowledge  about  a  source  vector  except  |or  the  i  nsh  mi  of  a  wavelront. 
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tested  by  the  fact  that  the  corresponding  ('SI)M  is  a  rank-one  mat rix.  This 
is  tin'  minimal  requirement  lor  spatial  processing  to  he  possible.  Intensive 
experimentation  at  sea  is  recommended  in  order  to  determine  the  extent  to 
which  this  last  assumption  can  he  trusted  across  various  aperture'  size's.  A 
signal  is  any  seuiree  with  a  wavefnuit.  a  sei-called  “perfect"  souree-.  Noise-  is 
everything  else.  We  are  doomed  t e »  simultaneously  eletermine  all  I  lie-  signal 
seuiree  vectors  in  a  sort  of  iouiijiiui  process. 

In  principle  the  elata  CSDM  measure-el  at  the-  sensor  outputs.  C, .  in¬ 
volves  several  "perfect'-  senirces  plus  noise.  The-  CSDM  e»f  P  nemceure-late-d 
perfect  source's  is  C,.  a  rank-/’  matrix.  If  P  we-re  large-r  than  .V  t  he-re- 
wenild  he  ue>  way  t < >  disentangle  the-  perfect -seuiree-  wave-fremt s  from  e-aeh 
other  (even  the  knowledge  e >f  C,  is  ne:>f  eneuigli).  For  a  solution  to  he-  pos- 
sihle  the  array  lias  t e >  he  large  enough  for  the-  data  CSDM  1  < *  he-  the-  sum  of 
thepe-rlect  seuirces  ( \S DM  and  e>f  t  lie-  noise-  CSDM ,  C,  C,  i  M.  whe-re-  t  lie 
tmise-  is  uncenre-late-el  with  any  perfect  seuiree-.  From  the  e-xpe-rime-ntal  elata 
ol  Cr.  the  first  problem  is  te»  de-rive  matrix  C,  and  therefore  the-  numhe-r  P 
e >1  perfee-f  source-s. 

!  his  coiilel  he-  eleuie-  accurately  at  the-  ce>st  eif  risky  he-tting  about  the- 
wavefreuits  i(  matrix  C,  is  known  to  iiave  a  specific  form.  Although  this 
is  sometime-s  ehnie-  in  the  literature,  we  have-  discarded  such  an  approach. 
Instead.  we-  prefer  Irving  to  eliminate  M.  Through  a  careful  examination  of 
the  eigenvalue's  of  matrix  C, ,  we  use-  experience  te>  cemject  ure-  how  many  of 
the-  largest  eigenvalue's  re-present  source-  ve-etors.  thus  determining  the-  value- 
of  /’.  I  lie1 2  me-theul  is  hacked  by  two  re-marks: 

I.  (blessing  abend  the  eigenvalues  is  less  risky  than  guessing  about  the- 
seuiree-  vecteirs;  since  every  semree  vecteir  elepends  on  all  the-  eligible- 
eigenvalues,  an  error  in  one  of  (hem  is  only  an  error  within  a  small 
part  of  the  final  source  vecteir. 


2.  Cuessing  is  le-ss  and  less  risky  when  the  array  greiws  in  size-. 

A  "large-,  ami  t  he-re-lore-  e-flieie-nl  at  l  av  has  a  wuh-  span  of  <-ige-n  value-s.  from 
the-  large-st  to  tin-  smalle-st.  he-|<uiging  to  the-  data  CSDM  with  I  lie-  smallest 
ih  in  asnuj  as  A  increase's.  I  bus  a  comparatively  vanishing  smallest  e-ige-n 
value-  is  a  sign  that  tin-  array  is  suf  lie  ie-nt  ly  large-  to  surmount  I  lie-  comph-xit  v 
ol  I  lie-  acoustic  f  i  e  •  |  <  I  wllii'll  IS  suppoSe-e||v  he-ing  eh-si'rihe-el  e  <  irre-et  I V.  I.xpe-ri 
men  I  s  at  se-a  a  re  r<-epiir<-d  l<  >t  <  •  I  i  « -  < '  k  i  n  c.  I  In-  <-x  I  «-nl  to  which  real  arra  vs  fulfill 
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this  condition  under  operational  conditions,  or  conversely,  what  size  array 
might  he  needed. 

II  tin'  span  of  the  eigenvalues  happens  to  he  only  a  guideline  for  the 
c  hoice  of  /’,  other  clues  have  to  he  used,  since  we  must  separate  matrix  Cr 
from  matrix  M.  As  described  here,  this  step  is  a  slightly  generalized  version 
of  some  methods  suggested  in  the  literature.  It  consists  in  representing  M 
with  a  reliable  general  model  including  free  parameters  which  are  adjusted 
to  bring  matrix  C,  to  rank  P.  If  the  noise  model  should  prove  to  be  highly 
reliable  (which  may  well  happen  for  fixed  arrays),  the  choice  of  P  may  be 
affected.  When  this  happens  the  free  parameter  adjustment  brings  matrix 
C,  to  the  minimum  possible  rank  compatible  with  the  noise  model. 

Besides  the  existence  of  the  wavefront  for  the  perfect  sources,  this  para¬ 
metric  noise  model  is  another  thing  which  must  be  restored  from  total 
ignorance.  Fortunately  the  final  choice  for  P  is  not  very  critical,  and  is 
subject  to  a  posteriori  verification  once  the  source  vectors  are  determined. 
Having  determined  the  number  ol  “  perfect"  sources,  F,  and  matrix  C*,  the 
next  problem  is  to  separate  the  wavefronts  from  each  other.  Apparently  it 
is  impossible  to  get  the  source  vectors  without  additional  assumptions,  and 
the  information  needed  corresponds  to  the  degrees  of  freedom  in  an  arbi¬ 
trary  (  P  I’)  unitary  matrix.  Any  larger  constraint  would  assign  matrix 
C,  specific  features  which  might  conflict  with  the  experimental  data  from 
which  C,  was  determined.  While  less  constraints  would  leave  the  problem 

unsolved. 

Actually  we  are  confronted  with  an  expression  for  the  source  vectors  in¬ 
volving  the  genuine  matrix  C,  (its  eigenvectors  and  eigenvlaues)  plus  some 
free  parameters  corresponding  to  the  degrees  of  freedom  mentioned  above 
in  (  1.21)  or  (  1.2o).  We  must  then  impose  constraints  which  are  just  fitted 
to  the  degrees  of  freedom  left  in  this  “CSDM-derived”  source  vector  expres¬ 
sion.  The  constraint  we  choose  requires  each  source  vector  to  be  expressed 
as  a  finite  coherent  sum  of  unknown  “conditional  plane  waves'’  CPW,  with 
1  he  number  ol  plane  waves  exactly  satisfying  the  available  degrees  of  free¬ 
dom.  So  we  get  a  “( '[’  W-derived  '  expression  for  the  source  vectors  in  such 
a  way  that  when  both  equations  are  satisfied  lor  the  whole  set  ol  source 
vectors,  we  find  a  sufficient  number  of  scalar  equations  to  compute  all  the 
free  scalars  in  both  expressions. 

The  ( '  1  *  V\  is  a  flexible  paramet  ric  ex pression  which  occasionally  proves 
to  In  a  plane  wave  across  the  uselul  bandwidth.  Although  difficult  to  carry 


out,  (his  way  of  determining  the  source  vectors  has  several  advantages: 

1.  It  preserves  the  experimental  data  and  depends  only  on  the  fitness 
of  a  parametric  noise  model  having  the  flexibility  usually  associated 
with  adaptivity. 

2.  It  only  imposes  constraints  on  the  source  vectors  which  are  exactly 
compatible  with  the  array  size  and  number  of  sources.  With  the 
source  vectors  at  hand,  an  approach  to  the  localization  problem  is 
briefly  outlined. 


Chapter  5 

GENERAL  CONCLUSION 
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In  chapter  2.  while  discussing  the  evolution  of  increasingly  complex  spatial 
processors,  we  introduced  the  concept  of  Directivity  and  its  natural  imple¬ 
ment  at  ion  via  classical  fixed  beam  formers.  This  classical  approach  requires 
two  stringent  concepts: 

1.  A  plane  wave  signal  coming  from  a  given  direction. 

2.  An  extremely  specific  and  unnatural  noise  model,  based  on  perfect 
isotropy  and  the  absence  of  correlation  between  any  two  nonoverlap¬ 
ping  elementary  solid  angles. 

I  he  above  leaves  an  insufficient  amount  of  information  to  exploit,  for  the 
"identiliation”  of  a  source  by  the  determination  of  its  wavefront  and  spectral 
density.  Here  the  wavefront  is  already  known  and  the  best,  signal  to  noise 
ratio  is  obtained  for  the  particular  sensor  combination  which  maximizes  the 
directivity  factor. 

When  trying  to  go  further,  through  an  estimation  of  the  “noise  alone” 
on  !!:<>  output  of  some  other  appropriate  sensor  combination,  we  fall  on  the 
rigid  concept  which  grants  no  specific  status  to  remote  sources  transmit¬ 
ting  plane  waves  lrom  directions  other  than  the  signal  direction;  remote 
sources  also  produce  noises,  or  “unwanted  signals"  on  the  “beam"  which 
we  would  like  to  clear  of  anything  but  the  signal.  Nevertheless  under  some 
assumptions  discussed  in  chapter  3,  it  is  possible  to  reach  the  notion  of 
“noise  alone"  and  its  estimation.  But  all  such  procedures  rely  on  an  excess 
ol  confidence  in  the  highly  artificial  notion  of  omnidirectional  noise.  In  fact 
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(lie  Directivity  concept  does  not  provide  reliable  removal  of  the  noise  spec¬ 
tra!  density  from  the  raw  data,  and  is  thus  unable  to  reach  the  realm  of 
the  so  called  “hyper-resolution"  systems  where  the  signal  spectral  density 
is  supposed  to  he  totally  rid  ol  noise. 

The  next  approach  to  spatial  processing  is  the  adaptive  beamiormer 
for  which  the  signal  wavelront  is  still  known  (its  spectral  density  needs 
to  he  identified)  but  lor  which  no  other  assumption  is  made  about  the 
unwanted  noise.  The  same  status  also  covers  sources  wit  h  a  wavelront  (but 
different  from  the  signal's)  and  diffuse  sources  with  no  wavefront;  it  is  all 
unwanted  noise,  (liven  the  signal  wavelront,  building  a  sensor  combination 
which  preserves  this  signal  (a  “signal  all  pass  filter ” )  is  easy,  but  the  output 
still  contains  some  unwanted  noise.  Knowledge  of  the  signal  wavelront 
also  makes  it  possible  to  build  sensor  combinations  which  do  not  let  the 
signal  through,  even  when  present-.  These  combinations  are  tin’  “Signal 
Tree  References**  (SKR).  Clenerally.  (A  I)  SFUs  are  available  from  an 
A -sensor  array,  and  there  are  several  ways  to  render  (hem  uncorrelated 
using  processes  which  are  the  physical  analog  of  the  diagonalizat ion  of  the 
SFR  cross-spectral  density  matrix  (C'SI)M).  Nevertheless  the  SKRs  remain 
correlated  with  the  noise  gathered  on  the  signal  output. 

In  this  paper  a  way  is  described  to  progressively  subtract  contributions 
from  the  noises  that  were  adaptively  derived  from  the  SKRs.  Strange  but 
t  rue.  I  hese  subt  ract  ions  are  hot  h  t  hose  of  noises  and  I  heir  spec!  ral  densit  ies. 
So  that,  in  a  repetitive  operation  exhaust ing  all  SKRs  the  residual  noise  left 
in  the  (unaltered)  signal  company,  has  been  brought  to  a  minimum.  This 
residual  noise  cannot  be  estimated  directly  since  we  know  neither  whether 
nor  when  a  signal  may  be  present.  Nevertheless  and  this  is  a  success 
of  the  adaptive  beamforming  theory  -it  is  possible  to  show  that  only  out 
scalar  relation  between  the  elements  of  the  CSDM  is  needed  to  obtain  an 
estimate  of  the  minimized  noise  above,  thus  reaching  the  realm  of  hvper- 
resolution  systems.  Such  a  relation  must  be  derived  from  intuition  or  more 
reliably  bound  to  a  particular  situation.  For  a  large  matrix  this  is  a  “weak" 
assumption  which  does  not  appear  to  be  loo  critical. 

Such  an  assumption  can  also  be  considered  as  automatically  pertaining 
to  "large  arrays.  '  In  chapters  !{  and  I  a  large  array  is  a  one  whose  size 
I  or  complexity,  or  number  A  I  is  large  enough  to  snrnioiinl  tin  coni j>h  .nl y 
of  the  number  of  sources  plus  the  complexity  of  the  propagation  medium. 
I  lie  test  of  this  capacity  is  derived  from  the  examination  ol  the  span  <>l 


( oxarsiox  73 

llic  (SIXM  eigenvalues.  in  every  frequency  bin.  It.  must  be  noticed  that, 
more  sensors  lowers  the  value  of  the  smallest  eigenvalue  and  is  liable  to 
widen  the  span  of  values.  It  is  important  for  the  smallest  eigenvalues  to  be 
much  smaller  than  the  largest.  We  then  can  consider  just,  neglecting  this 
smallest  eigenvalue,  which  immediately  provides  the  wanted  scalar  relation 
mentioned  above  and  yields  the  signal,  in  principle,  rid  of  noise  inasmuch  as 
second  order  statistics  are  concerned.  This  “large  array  condition”  remains 
to  be  experimentally  checked  on  actual  arrays  in  operational  conditions. 
Such  a  check  must  done  before  stating  that  the  array  correctly  represents 
the  medium  and  sources,  and  therefore  is  an  efficient  array. 

Noise  subtraction  techniques  and  noise  cancelling  techniques  are  all  par- 
t  icular  branches  of  t  his  present  at  ion  as  is  the  removal  of  self-noise  from  hull- 
mounted  arrays.  Although  one  may  be  enthusiastic  about  the  idea  of  nearly 
achieving  hyper-resolution  while  dropping  all  assumptions  concerning  un¬ 
wanted  noise,  some  caution  is  in  order  because  the  process  is  very  sensitive 
to  the  shape  of  the  signal  wavefront.  A  small  deviation  of  a  signal  away 
from  its  expected  wavefront  unpredictably  alters  the  output  signal  spectral 
density,  a  quantity  we  had  been  careful  to  preserve.  Making  the  process 
more  robust  is  possible  but  this  incurs  a  loss  in  the  noise  reduction  capac¬ 
ity.  Besides  it  is  impossible  to  be  prepared  for  all  the  possible  wavefronts 
describable  in  an  .Y-dimensional  space  (a  generalization  of  beamforming), 
io  do  so  would  be  both  fantastically  complex  and  a  true  waste,  since  we 
are  i >nly  considering  “large  arrays,"  with  far  fewer  sources  than  sensors. 

Finding  a  way  to  process  signals  with  unknown  wavefronts  forces  us 
to  restore  some  “weak"  assumption  the  weakest  possible.  We  base  our 
a -m ti i pt  ion  on  the  observation  that  a  source  with  an  unknown  wavefront 
iimo/  m  n  rlln  It  ss  hurt  somt  “existing"  wovefront.  Whether  this  condition 
is  fulfilled  in  the  real  world  will  be  decided  by  experimental  studies  which 
determine  whether  an  “interesting”  source  of  the  type  we  might  have  to 
identify  exhibits  a  rank-one  CSDM  regardless  of  the  wavefront.  This  is, 
ah  >ng  wit  h  I  he  eigenvalue  span,  t  he  corner-stone  of  spatial  processing  and  its 
ultimate  justification.  One  may  expect  the  corner-stone  is  firmly  attached 
to  reality.  especially  lor  large  fixed  arrays. 

<  ) 1 1  the  other  hand,  the  global  ('SUM  of  several  "perfect."  sources,  each 
po-sessing  a  wavefront .  is  not  by  itself  sufficient  for  t  he  separat  ion  ol  several 
M.nire  vectors.  ( )  1 1  e  c  a  1 1  on  ly  sa  v  /  Ik  run  I'  of  1 1)  is  “pt  r  ft  el  so  u  ret  s  CSI)M"  is 
I  In  n  ,t  in  In  r  I '  of  sou  ret  s.  f  u  1 1  lieniiore  noise  sou  ices  wit  li  a  dill  use  st  met  ure 
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(not  granted  a  wavefront)  add  their  own  full  rank  ('Sl)M  to  yield  t lie  only 
accessible  matrix,  the  raw  data  CSDM.  The  procedure  to  he  considered  is 
t  hen  a  t  wo-st  ep  one. 

The  first  step  consists  in  trying  to  eliminate,  as  much  as  possible,  the 
noise  ('SI)M  from  the  raw  data  CSDM.  To  do  so  we  need  a  reasonable 
general  model  o!  Ine  noise  CSDM.  which  is  made  flexible  enough  through 
the  use  of  free  parameters  (in  principle  (A  /’)  parameters  if  /’  is  already 
chosen),  to  be  fitted  to  the  widest  category  of  situations.  Then  the  free 
parameters  are  varied  to  reduce  the  rank  of  the  difference  between  the  raw 
data  CSDM  and  the  noise  CSDM  from  A  to  /’.  The  remaining  difference 
matrix  is  expected  to  be  the  source  CSDM. 

While  the  number  of  /’  sources  might  usually  be  derived  from  an  ex¬ 
amination  of  the  sef  of  eigenvalues  in  the  raw  data  CSDM.  under  some 
conditions  the  noise  model  may  be  reliable  enough  (especially  for  fixed  ar¬ 
rays)  to  either  impose,  or  at  least  influence  the  choice  of  /’.  In  fact  the 
choice  ol  /’  cannot  be  strictly  determined  precisely  because  we  have  made 
no  assumption  about  the  wavefronts.  On  the  other  hand,  this  is  not  a 
critical  choice  since  it  is  subject  ♦<>  a  posteriori  verification. 

In  the  second  step,  /’  is  now  known  and  the  second  order  statistics  of 
the  noise  have  been  eliminated.  What  remains  is  to  disentangle  the  source 
vectors  from  their  (A  •  /’)  CSDM.  Unfortunately  this  disentanglement 
cannot  be  don*'  without  some  additional  assumptions  because  there  are 
an  unlimited  number  of  solutions  lor  source  vectors  directly  derived  from 
the  CSDM.  These  solutions  depend  on  a  set  ol  undetermined  parameters 
corresponding  to  the  degree  of  freedom  ol  a  ( /’  •  /’)  unitary  matrix.  As 
a  consequence  we  have  to  impose  some  "weak"  constraint  on  the  source 
vectors,  so  that  the  lree  parameters  become  computable  thus  allowing  a 
solid  ion  to  emerge. 

I  he  constraint  suggested  here  is  a  parametric  representation  of  a  source 
vector  as  a  coherent  sum  of  a  limited  number  of  “conditional  plane  waves," 
CDW,  which  are  not  exactly  plane,  but  whose  geometry  is  only  subject 
to  a  posteriori  control  while  the  parametric  value  remains  constant.  The 
limited  number  ol  CPW  per  source  vector  ami  their  flexibility  in  terms  of 
physics,  seem  to  justify  the  claim  lor  the  assumption  being  “weak." 

I  lie  constraints  on  the  CPW  in  this  source  vector  model  are  chosen  in 
such  a  way  that  the  number  of  free  scalars  on  both  sides  exactly  match 
the  number  of  relat  ions  old  aim'd  wlun  I  lie  "CSD.M  derived"  expression  is 
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confronted  with  the  same  source  vector.  Most  of  the  assumptions  about 
the  source  vectors  which  have  been  proposed  in  the  literature  are  far  from 
this  point  of  balance,  and  actually  impose  a  specific  form  onto  an  already 
known  experimental  C'SDM. 

Our  notion  of  CPW  seems  to  be  acceptable  as  a  weak  assumption  be¬ 
cause  in  a  rich  parametric  representation,  the  physical  definition  of  a  given 
parameter  tends  to  be  blurry  without  loss  ol  the  parametric  value.  Also, 
this  might  be  true  for  a  parametric  localization  process  for  the  assessment, 
of  available  free  parameters  in  the  description  of  the  medium,  such  as  the 
one  briefly  sketched  at  the  end  of  chapter  1. 
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Appendix  A 

Classical  Relation  between 
Cross-Spectral  Densities 
Before  and  After  Filtering 


flic  following  well  known  rule  is  one  of  (lie  basic  tools  of  flic  mathematical 
present  at  ion : 

li\  -i<l>,i . D* 


IK  -|<Kj  •-  B" 

(liven  two  input  noises  H j  am!  U i  respectively  filtered  with  frequency 
responses  'l> i  and  'In  which  deliver  two  output  noises  /l,’  and  IK.  we  have 


csd  [u\lk]  v 


(A.I) 


where  csd  denotes  cross  spectral  density  and  sd  spectral  density.  In  the 
particular  case  where  U\  and  U\  are  the  same  noise  then. 


csd  [/*"/?“]  sd[R I't’pl', 


•WWW 


SW! 


V  ^  ^ 


Appendix  B 


Matrix  of  Filters  and  Related 
Equations 

As  ;i  rehersal  with  the  CSDM  aiul  a  matrix  of  filters  let  us  consider  N 
original  noises: 

/?!.  B2,  .  ■  ■  Bn - Bm , . . .  Z?/v : 

with  a  CSDM  C  whose  elements  are  c„„,  (;?th  row,  mUl  column).  Let  these 
elements  he  processed  through  a  matrix  of  filters  as  shown  in  figure  13.1. 
The  frequency  response  of  a  filter  is  gpn  where: 

I.  I  he  first  subscript  refers  to  the  p'h  output  noise  Bp  which  is  filtered 
by  t  he  />,h  row  of  filters. 

The  second  subscript  refers  to  the  u"’  input  noise  B„  which  feeds  all 
the  filters  of  the  n,h  column  of  filters. 

So  the  matrix  of  filters  G  appears  as  a  matrix  is  usually  written — with 
elements  in  the  plU  row  and  the  77 th  column.  Using  the  symbol  <g>  for 
“filtered  by",  we  have: 

It"  i  B,  //;,i )  -1  ...  1  I  B„  (?:  g,„, )  T  •  •  •  T  (  B,„  0  gpm  )•!...  +  (  By  CO  gpN ) , 

IB  ( /f,  //,,,  )!.■•*(  B„  W  )*-...  1  (  B,„  to  )  t  . . .  4  (  By  0  <i,hy  ) . 

Using  the  rule  from  appendix  A.  the  cross-spectral  density  (csd)  between 
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INPUT  NOISES  (CSDM  :  C) 


\  A' 

!lp<)  ^  ^  (B.l  ) 

»,  1  hi  1 

which  is  ail  element  of  the  output  ('SI)M  Y. 

Now  q'  is  the  element  of  the  matrix  G  ,  so  that  i/m  can  he  written 

N  .V 

ji  yn]  ^  ^  imil  nit/  (  H .  t 

n  I  ■-.  I 

which  is  exactly  the  mathematical  expression  o|  I  lie  element  for  a  product 
of  three  matrices.  Indeed  an  element  of  H  GCD  is 


t  BP  END  LX  B 


XI 


Comparing  (133)  with  (132)  clearly  shows 


Y  -  GCGll 

i  he  spectral  density  of  /C  is 


N  N 

^  ,  9p» r"  mfjpm- 


(13.4) 


(13.5) 


Considering  the  sequence  of  gpn  with  fixed  />  as  the  row  vector  g;),  the 
sequence  of  <y'mi  becomes  t  he  column  vector  g|.  Therefore 


y„  =  gpCgt.  (B.6) 

the  quadratic  form  of  vector  g;)  (row  matrix)  through  matrix  C. 

To  follow  the  convention  made  in  section  3.3,  we  should  represent  the 

set  ol  litters  g;)  by  a  column  vector  V  —  g|  so  that  t ipp  —  V^CV  as  in 
section  3.1  equation  (3.13).  and  section  3.5. 
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Classical  Factorization  of  a 
CSDM  -  Normalized 
Eigenvectors 

A  positive  definite,  llerniitian  matrix  like  the  CSDM  C  is  known  to  factor 
into  t  he  product 

C  -  UEljt.  (CM) 

where: 

E  is  t  he  diagonal  mat  rix  of  t  lie  cigenvlues,  <  i .  e2, . . .  < 

U  is  the  unitary  matrix  whose  column s  are  the  normalized  eigenvectors  of 
C:  i i i .  lij. .  .  .  u,v;  placed  in  the  same  order  as  the  eigenvalues  in  E. 

Ih'cause  these  eigenvectors  are  normalized  and  orthogonal  to  each  other: 


5: 

V. 


■>. 


li 0.  (/>  A  </) 

UTU  Iv.  the  identity  matrix. 

One  may  rapidly  check  that  equation  ((’.It  is  the  right  one  when  the' 
columns  ol  U  are  the  eigenvectors.  The  definition  ol  li,,  is. 
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If  equation  (CM)  is  true  we  must  have 

UEUtup  -  tpuP-  (0.2) 

We  cheek  this  by  multiplying  out  the  left  hand  side.  Beginuing  with  U^ii,,. 

and  noticing  the  rows  of  U  are  lip,  we  end  up  with  a  column  vector  v 
where  all  the  elements  are  null  except  the  p' '  whose  value  is  one.  Pursuing 
with  Ev,  we  get  a  column  vector  with  all  elements  null  except  the  ptU  whose 
value  is  fp.  Finally  in  the  matrix  product 


only  the  elements  of  the  plh  column  of  U  are  scaled  by  the  common  factor 
i  r.  The  result  is  then  (pup  as  expected  from  ((’.2). 


Appendix  D 

Orthogonal  Images  of  a  Set  of 
Correlated  Noises 


We  now  know  that  I  h rough  t  he  matrix  of  filters  G,  I  lie  former  CSDM  C 
has  liecome  GCG^  which  may  also  be  written  as  GUEU^G^.  1 1  we  choose' 
to  use  a  matrix  of  filters  G  -  U^,  clearly  the  out  put  matrix  becomes 


U^ul  E  IVu 


(D.l) 


lnrausc  U'U  I  v  is  the  identity  matrix.  So  we  have  got  a  set  of  nncor- 
related  noises  with  spectral  densities  1 1 ,  c2. . . .  f)M  .  . .  t:) . . .  <■  ,\- .  'This  is  the 
physical  interpretation  of  the  diagonalizatiou  of  matrix  C.  'I'lie  new  inde¬ 
pendent  noises  are  the  '‘orthogonal  images"  of  the  input  noises.  The  choice- 

ol  G  U  1  shows  that  a  row  of  filters  in  G,  is  the  corresponding  row  of 

the  Ut  elements,  and  therefore  the  complex  conjugates  ol  the  elements  of 
the  corresponding  column  in  U.  In  oilier  words  the  normalized  eigenvector 
U.  is  the  column  vector  which  is  representative  of  the  set  of  filters  made 
of  the  complex  conjugates  of  its  elements.  We  are  therefore  faithful  to  the 
convention  staled  in  section 
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Appendix  E 

Canonical  Form  of  a  CSDM 


I  here  is  mint  her  canonical  lonn  of  matrix  C  which  is  of  use  in  chapter  ■!. 


C  V  £„  U„tljl  . 


(E.l) 


In  I  (rackets  is  the  dyadic  product  of  a  normalized  eigenvector  which  forms 
an  (A  •  A  )  matrix.  It  is  easy  to  check  the  previous  equation,  for  example 
with  111  e  eigenvector  li t , 


Cii,  (»!»!  )  < 

»  -  1 


■re  the  scalars  are  between  parenthesis.  Since 


ii'Ci]  -  1.  and 


■Ve  end  up  wit  I) 


ii,,Ui  -  0.  for  n  ~  1, 


C  fi  i  ([ii|. 


which  is  the  very  definition  <>(  t  i  ( .  I  he  same  then  holds  true  for  all  ii(>. 
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Section  4.6 

1  cl  us  consul 

er  P  vectors:  V|. 

...Vr.  with  /’  A'  des. 

•rihiiu;  a 

true  /  ’ 

dimensional  s 

tihspace  (  no  linear  i 

dependence  het  ween  any  | 
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.in  A -dimensional  vector  space. 

Also  consider  a  diagonal 

(A  .V) 

mat  rix 

II  [see  equation  (-1.1  1  )|.  where  the  first  P  diagonal  elements  are  1. 

and  the 

remaining  (  A 
Ui. 

/’ )  others  are  0 
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equation  f  1. 12)  l>y  V [ , 


r  , 

iiV,  -  v  v,,  (vjv,) . 

1 

In  this  relation  V,  appears  to  he  a  linear  comhinal ion  of  the  other  vectors 
V(>.  Hut  this  is  impossible  since  they  have-  been  taken  as  describing  a  (me 
/'-dimensional  snbspaee.  The  only  wav  to  comply  with  this  equation  is  to 
accept : 

V^V  ]  -  0  for  p  /■  I . 

vfv,  -  I. 

W'ctor  V[  is  a  normalized  vector  orthogonal  to  the  others.  YVhal  is  true  for 
V] .  is  true  for  any  V,,.  Therefore  feature  2:  the  V()  are  normalized  vectors 
ort hogonal  to  each  ot her. 
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